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For centuries, military forces have used camouflage to obscure potential targets 
from the enemy. Because the eye is fairly adept at picking out edges, colors, and 
bright areas, camouflage is often used to degrade these qualities from human detection. 
The purpose of this thesis was to investigate the role of certain spatial, temporal, and 
chromatic features on the human visual system and how these features may aid the 
quest for better camouflage. Methods: Test patterns were spatio-temporal raised 
cosines of varying orientation (horizontal or vertical and oblique), spatial frequency (1, 
3, and 7 cpd), and modulated at 2.0 Hz.- Color contrast thresholds were determined 
from 16 different red-green color mixture ratios. This methodology eliminates the 
problems with luminance artifacts and the need to determine exact equiluminance. 
Results: The data formed an ellipse with the half-length measuring color discrimination 
and the half-width measuring brightness discrimination. A maximum likelihood 
method was used to fit the data. Three of the four subjects showed a 3 cpd chromatic 
oblique effect, while the 1 and 7 cpd achromatic and chromatic oblique effect was 
inconsistent across subjects. Conclusions: While real-world objects are more complex 
than laboratory stimuli, knowledge of spatial and chromatic qualities that inhibit 


detection will aid the quest for better camouflage. 
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EXECUTIVE SUMMARY 


For centuries, military forces have used camouflage to obscure potential targets from 
the enemy. Because the eye 1s fairly adept at picking out edges, colors, and bright areas, these 
are the qualities that camouflage seeks to degrade. The purpose of this thesis was to 
investigate the effects of certain spatial, temporal, and chromatic features on the human visual 
system and how these features may aid the quest for better camouflage. 

Camouflage targets tend to match their backgrounds both in: color and structure. 
Netting or paint can make detection of a potential target more difficult by reducing chromatic 
contrast. They can also reduce structural contrast by reducing sharp edge effects. This is 
evident in the "stealth" design in which the lack of defined edges helps reduce the radar 
signature. By removing defined edges, and thus high spatial frequency components, the visual 
signature is also reduced since "stealth" ships and aircraft generally have backgrounds 
consisting primarily of low frequencies (sea and sky). The perception aspects of this real world 
example can be simplified by examining less complex stimuli in the laboratory. For example, 
camouflage design would be affected if it were known that oblique chromatic lines were more 
difficult to detect than honzontal or vertical chromatic lines. 

It has been shown that honzontal and vertical achromatic lines are easier to detect than 
oblique achromatic lines. This phenomenon is known as the "oblique effect," which in this case 
is an achromatic oblique effect. Several studies have shown that the magnitude of the oblique 
effect is largest with high spatial and a low temporal frequency sinusoidal gratings. Previous 
researchers have used this knowledge to design experiments testing for a chromatic oblique 
effect, but have had problems with luminance artifacts due to the difficulty of obtaining exact 
equiluminance. By adapting the methodology of an earlier researcher the problem of 
determining exact equiluminance was avoided and an experiment to test for a chromatic 
oblique effect was designed. 

The experiment was conducted concurrently at the Naval Postgraduate School (NPS) 
and the University of Louisville, Kentucky (UL). Four subjects, 2 NPS and 2 UL, volunteered 


for this experiment. All subjects had normal (20/20), or corrected to normal, acuity and color 
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vision. Stimuli were presented by a VisionWorks computer graphics system (Vision Research 
Graphics, Inc.) on an IDEK MF-8521 high resolution color monitor (21" X 20" of viewable 
area). The monitor had a resolution of 800 by 600 pixels (x=75.02 and y=74.92 pixels/degree), 
98.9 Hz frame-rate, mean chromaticity of r = 0.334, g = 0.336, b = 0.300 (1931 CIE), anda 
maximum luminance of 100 cd/m2. The University of Louisville's apparatus and procedure 
were identical to the Naval Postgraduate School's, except that the stumuli were displayed on a 
17" Nanao Flexscan F2.21 color monitor. Subjects viewed the monitor from 1.5 meters and 
were positioned by an adjustable chinrest. 

Sinusoidal gratings were presented within a spatially windowed circular test field that 
subtended 7.59° of visual angle. The Gaussian window was truncated at +1 standard 
deviations for both x and y directions. The test patterns were one-dimensional spatio-temporal 
sinusoids of varying orientation (principal and oblique), spatial frequency (1.0, 3.0, and 7.0 
cycles/degree), and color contrast. Test patterns for each subject consisted of two onientations, 
principal (0° and 90°) and oblique (45° and 135°). For each subject, maximum sensitivity for 
each onentation within the principal and oblique grouping was chosen. All sinusoids were 
raised cosines temporally modulated at 2.0 Hz. The sinusoid pattem was presented in a 1500 
msec interval with contrast ramped on and off according to a linear window. (Contrast 
peaked at 202 msec and fell at 1304 msec. Color contrast was computed by different ratios of 
percent red and green. Sixteen different sinusoidal red-green color mixtures were generated by 
changing the red phosphor only, green phosphor only, or by changing the red and green guns in 
fixed proportions. Color contrast was defined according to the Michelson formula. Thresholds 
were determined by a sequential two-alternative forced choice adaptive psychometric 
procedure, QUEST. Threshold was defined at 75 percent correct. A total of 480 tnals, 30 
tnals per condition, were randomly presented within each session. A session (~ 45 minutes) 
consisted of one sinusoidal condition with 16 different red-green color mixtures. A subject had 
to complete six sessions to contribute one set of 16 thresholds for each condition. 

Numerous surveys of differential thresholds have been carried out, but one of the more 


extensive ones was completed by D.L. MacAdam. The data from this survey was elliptical (the 


closed curves connecting the thresholds were elliptical in shape). It was shown that the errors 
of thresholds about these closed curves were Normally distributed, therefore the curves should 
be ellipses, as they appeared to be. 

The elliptical properties of the experimental data were used to fit ellipses using the 
method of maximum likelihood a nested F-test (which is only approximate, because of non- 
linearity) was used to examine the significance of the different orientations. Results of this test 
showed that the chromatic oblique effect was inconsistent across the spatial frequencies of 1, 3 
and 7 cycles/degree (cpd). The lack of a 1 cpd chromatic oblique effect was due to the 
insensitivity of the chromatic channel at the lower spatial frequencies. Two of the four subjects 
showed a graphical 7 cpd chromatic oblique effect, but this was non-significant. Chromatic 
aberration may have been a factor. Three of the four subjects showed a 3 cpd chromatic 
oblique effect with two significant and the third marginally significant. It is predicted that the 
marginally significant subject would show significance with additional trials. 

The main value of this study 1s the tool tt provides for further investigation of a 
chromatic oblique effect without the problems associated with luminance artifacts. 
Additionally, further investigation of a chromatic oblique effect will likely provide knowledge 
of spatial and chromatic qualities that inhibit detection. Knowledge of these qualities will aid 


the quest for better camouflage 
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I. INTRODUCTION 


The continual decline of the military budget has necessitated the increased protection of 
current and future war-fighting assets. This increase, coupled with public expectation of zero 
or close-to-zero casualties, has forced the services to reassess the way they conduct operations. 
Today's political climate necessitates this reassessment since the potential loss of public support 
for military actions generally increases with the number of American casualties. 

A cruise missile attack is an effective method used to minimize civilian and friendly 
casualties. For example, cruise missile attacks were successful against Iraq’s military targets 
during the Gulf War. Although the cruise missile is an effective weapon, it is an expensive 
resource for the United States military inventory. A less expensive alternative is an air strike, 
but the disadvantage of an air strike is the increased likelihood of aircrew casualties and missed 
targets. Ideally, the aviator would like to enter the threat zone undetected, thereby increasing 
the probability of locating the target without becoming engaged by the enemy. 

The ability to avoid detection is a distinct advantage during battle. For centuries, the 
element of surprise has resulted in a quick and decisive destruction of forces. For example, the 
U.S. Air Force F-117 “stealth” fighter was responsible for much of the precision bombing 
during the Gulf War. Because their aircraft was nearly invisible within the Iraqi air defenses, 
pilots had additional time to accurately drop bombs on target. This near-invisibility must 
extend beyond the visible spectrum since today’s battlefields are equipped with electro-optical 
sensors that often extend the range and increase the probability of detecting potential threats. 
Electro-optical sensors make detection possible through visual, infrared, thermal and other 
means. For example, shielding hot parts of a vehicle is part of the vehicle's thermal camouflage 
as it tries to disperse its thermal signature. 

But, even with these technological advances, a large threat to operations, and to 
reconnaissance operations in particular, is often the most common sensor--an enemy’s eyes. 
Prior to the start of the Gulf ground war, U. S. forces had reconnaissance teams operating in 
the interior of Kuwait. These teams used their speed and camouflage to prevent detection and 


capture. If their camouflage had been inadequate, capture and death would have been likely, 
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delaying the start of the ground offensive. While the saying “if they can't see you, they can't 
hurt you” is no longer true, the likelihood of being hurt is reduced when the enemy cannot see 
you. 

To make the enemy’s task of detection more difficult, camouflage is used. Since 
detection has traditionally been associated with visual detection, camouflage is generally 
thought of as making the visual detection of personnel or any potential target more difficult. It 
is in this sense that camouflage will be discussed. 

Camouflage has many different applications, ranging from the clothes and face paint of 
an infantryman to the netting used to cover tanks and vehicles to the paint used on larger 
platforms such as aircraft and ships. Since the eye is fairly adept at picking out edges and 
bright areas, camouflage is often used to break up edges and to cover or conceal bright areas. 
All objects possess certain unique qualities of shape and color. In order to deceive a sensor, 
the object must blend with the background. An invisible object would match the background, 
while an easily identifiable target would contain noticeable spatial, temporal, or chromatic 
features. By manipulating the spatial, temporal and chromatic features, an object can be made 
more difficult to detect. Knowledge of the range of these features will aid military designers in 
their quest for better camouflage. 

To understand how to make detection more difficult, one must also understand the 
sensors that will be used. Since, in this discussion, the primary sensors are the enemies’ eyes, 
either directly or through some sort of image intensifying mechanism, a working knowledge of 


how the eye works and what cues it uses to accomplish detection is essential. 


Il. BACKGROUND 


The ability to perceive an object is an unconscious, automatic process. The world is 
filled with a variety of sensory information that stimulates our senses. This information is 
obtained through visual, auditory, tactile and olfactory inputs. We use our senses to collect this 
information and translate it into meaningful units of sensory awareness. This information is 
then relayed to the brain, resulting in the formation of perceptions. The brain then categorizes 
this sensory data and compares it to past experiences. Thus, perceptions are a culmination of 
sensory inputs that are organized into a meaningful representation of the outside world. The 
study of perception involves a complete understanding of the description of objects, 
appearances and events in the outside world (Sekuler and Blake, 1994). In brief, sensation and 
perception refer to a sequence of events: stimulation of an external object; machinery to 
capture this information; and translation of this information into electrical energy to form an 
experience. Perceptual experiences guide our actions in the world. This thesis investigates the 
different techniques of manipulating visual sensory input to change the appearance or 
perception of an object. (Sekuler and Blake, 1994) 

The basic function of the eye is to capture visual sensory input or light and focus it on 
the retina, a thin layer of receptor cells located in the back of the eye. The light must pass 
through these retinal cells before reaching the photoreceptors, which convert the light into an 
electrical signal. This process of converting an external stimulus into a neural signal is called 
transduction. These neural signals are then sent through the network of retinal cells, which, in 
tum, are sent to the brain (Figure 2.1). 

The human retina contains two major classes of photoreceptors, rods and cones. There 
are approximately 8 million cones and 120 million rods. The fovea, located at the center of the 
retina, has the greatest concentration of cones. The cones are sensitive to daylight and provide 
high-acuity color vision, while the rods are used for mght viewing and are thought to be 


achromatic. Between five and ten percent of the total cone population are sensitive to short 
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wavelengths, with a peak sensitivity of 420nm. This cone class is referred to as S or blue 
cones. A few S cones are located at the fovea, but the largest concentration of S cones forms a 
ring around the fovea. This large concentration of cones tapers off with increasing distance 
(eccentricity) from the fovea. The remaining cone population is sensitive to middle (peak 
sensitivity of 530nm) and long wavelengths (peak sensitivity of 565nm). The long-wavelength 
(L or red) cones outnumber the middle-wavelength (M or green) cones two to one. The R and 
G cone distribution is randomly mixed throughout the retina, with the greatest concentration 
located inside the fovea. (Tovee, 1996) 

While he did not directly identify cones, Thomas Young believed that the retina 
contained three receivers that were sensitive to a limited number of light vibrations. These 
receivers are now known as cones. Young also is responsible for one of the first explanations 
of how we perceive color. Earlier, Newton had demonstrated that white light could be split by 
a prism into a spectrum of colored lights. He found that recombining some of these colored 
lights resulted in the onginal white light. Newton mixed and subtracted colors, but he did 
not attempt to explain how we perceive color. Young, however, postulated that color 
perception is due to the vibrations of light interacting with the retina. The three receivers in the 
retina were broadly tuned with overlapping sensitivities. Helmholtz confirmed and elaborated 
Young’s color theory by showing that there are three types of receivers (cones or 
photoreceptors) in the human retina, and that each type contains a different pigment. The 
spectral sensitivity of the cone is determined by the absorption spectrum of its photopigment. 
(MclIlwain, 1996) 

Young hypothesized that there were three broadly tuned receivers in the retina 
because, he reasoned, a single broadly-tuned receiver could not provide enough information 
about the wavelength of light. If there were two receivers, then there would be one particular 
frequency that excited both receivers equally, and thus white light would be produced 
(intersection of the curves in Figure 2.2b). Since Young did not observe white light in the 
color spectrum produced by a prism, he concluded that the visual system must have three 
broadly tuned receivers. A single cone pigment cannot discriminate between changes in 


wavelength and changes in the intensity of light. A cone can only increase or decrease its 


output, so its signal is ambiguous as to whether the change is due to a shift in the light’s 
wavelength or to a change in its intensity. This type of response is explained by the principle of 
univariance. For example, a retina that contains one cone class may give the same response to 
two different wavelengths, while a two-cone class retina may excite the receptors in different 
ratios (Figure 2.2). Therefore, the retina needs at least two cone types to distinguish between 
changes in wavelengths. Pnmate vision can discriminate millions of colors with three cone 
types, whereas non-primate mammals, which usually have two cone types, cannot. These 
mammals tend to rely much more heavily upon auditory and olfactory sensory inputs to 


survive. (Tovee, 1996) 
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Figure 2.2. Wavelength discrimination by (a) a one-class retina and (b) a two-class retina. 
A retina that contains one cone pigment responds more or less with the same energy for 
some wavelength within its spectral sensitivity--the principle of univariance. A retina that 
contains two cone pigments will have different responses depending upon the location of 
the two wavelengths. From Mcllwain [ 1996]. 


The Young-Helmholtz tnchromatic theory accounts for many, but not all, of the 
phenomena associated with color vision. The theory predicts that a signal comprised of a 
certain combination of long and medium wavelengths cannot be distinguished from a specific 
third wavelength (yellow), but it does not account for the fact that the signal appears yellow 
rather then red-green or green-red. Additionally, the phenomenon where an object’s color 
appears to vary depending on the colors viewed immediately before viewing the object 
(successive color contrast) or on the colors surrounding the object (simultaneous color 
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contrast), cannot be accounted for. A theory that does account for both these phenomena and 
those explained by the Young-Helmoltz theory is known as the opponent-color theory. 
(MclIlwain, 1996) 

The opponent-color theory was first introduced in 1878 by Hering and has been 
furthered through the work of Hurvitch and Jameson. Hering postulated that there were three 
visual processes, two chromatic and one achromatic. The processes consist of three 
antagonistic or opponent pairings. These pairings are red-green and yellow-blue for the 
chromatic processes and black-white for the achromatic process. Such opponent pairs are well 
explained by the center-surround or on-center and off-center type ganglion cells consisting of 
the aforementioned pairings. These opponent pairs account for the fact that the colors in these 
pairings cannot be seen at the same time; thus, there are no reddish-green hues. The inputs 
from the S, M and L cones in the first diagram of Figure 2.3 display in a simplified way how 
these inputs are combined to form a signal that we perceive as blue. The second diagram in the 
Figure is identical to the first, except that the weightings of these signals result in the perception 
of yellow instead of blue. Intermediate weightings of the signals displayed result in our 


perception of many more colors than the six shown in Figure 2.3. (Tovee, 1996) 
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Figure 2.3. Simplified hypothetical display for a model based on color opponency. From 


Mcllwain [1996]. 


The cones relay their information by synapsing onto ganglion cells, whose axons travel 
through the optic nerves to the visual cortex located in the back of the brain. The human eye 
contains approximately one million ganglion cells. The input of 128 million photoreceptor 


signals has been reduced to an output produced by these one million ganglion cells. In 1938, 


7 


Hartline discovered that the retinal ganglion cells of the frog were compnsed of two 
concentrically shaped ring-like areas on the retina (Figure 2.4). He found that certain ganglion 
cells increased their electrical energy when a light was passed through their center or inner ring 
and decreased their electrical energy when a light was passed through the outer rng. Hartline 
called these ON-center ganglion cells. Cells that respond vigorously to a dark center and a 
light edge are called OFF-center ganglion cells. A network of such ON-center and OFF-center 
cells is responsible for providing edge detection as well as orientation information. Hartline’s 
work laid the foundation for later work by Devalois (1958), who discovered that primate 
ganglion cells behaved in a similar manner. (Sekuler and Blake, 1996) 

The morphological and physiological properties of primate ganglion cells are 
divided into three categories: large size P, or A cells, small size Pg or B cells, and P, or W- 
like cells. Livingstone and Hubel (1988) classify two major types of ganglion cells, P, and 
Ps cells, and their projections to different cortical regions within the primate visual system. 
P, cells (ten percent of ganglion cells) have high conduction velocities, transient 
responses, no color sensitivities, high contrast sensitivity and very good temporal 
frequency modulation. Pg cells (80 percent of ganglion cells) have lower conduction 
velocities, sustained responses, low contrast sensitivity, color opponency and moderate 
temporal resolution. Both P, and Pg neurons are segregated into two different pathways 
which project to different locations within the lateral geniculate nucleus (LGN), V1, and 
higher cortical regions within the primate cortex (Livingstone and Hubel, 1984). Refer to 
Figure 2.5 for a graphical representation of the visual pathway. (Merigan, 1989) 

The magnocellular pathway (M pathway) receives input from P,, ganglion cells that 
project first to layers 1-2 of the LGN. Cells in the magnocellular geniculate layers project 
to layer 4C, of the primary visual cortex. From layer 4C, they project to layer 4B, which, 
in turn, projects to visual area 2 and to the medial superior temporal area (Tovee, 1996). 
This pathway is thought to involve spatial awareness, that is, ‘where’ an object is located 
in space. Alternatively, the parvocellular pathway (P-pathway) receives input from Px 


ganglion cells and first projects to layers 3-6 of the LGN. Cells in the parvocellular 





at aa a etal 








ee 
208. 


{x 


“ta 


. 


| ganglion cell with (a) an “ON” center 


ina 
ibitory and (b) an “OFF” center excitory and “ON” 





£ 
Vv 


as 


$ 


Receptive field layout of a ret 


d “OFF” periphery 1 
periphery inhibitory. From Schiffman [1996]. 


of 


L. 


Figure 2.4 
excitory an 





Washington 


? 


laboratory home page 


S 


d from Dr. Van Essen’ 


aine 


Figure 2.5. Images were obt 
University, St. Louis, Mo. 


geniculate layers project to layer 4Cgin the primary visual cortex. From layer 4Cg they, 
in turn, project to layers 2 and 3. From these two layers, information is sent to visual area 
2, layer 4, and then to the inferior temporal cortex (Tovee, 1996). The inferior temporal 
cortex is thought to be concerned about the ‘what’ of an object. In summary, a crude but 
simple classification of the P and M pathways can be characterized as the ‘what and 
where’ of objects that an observer perceives. (Livingstone and Hubel, 1988) 

Scientists have leammed about the P and M pathways and their contributions to vision 
through studies that observed one pathway after the other pathway had been made inoperative 
by lesioning it. In 1990, Schiller and Logothetis created lesions in the P or M pathway of 
monkeys and then conducted various tests, including color perception, texture perception, 
acuity, pattern perception, flicker perception, and contrast perception. The results of these 
tests indicate that particular functions do tend to be associated with a specific pathway. In 
general, the P pathway is associated with color perception, texture perception, pattern 
perception, acuity, and contrast perception, whereas the M pathway is associated with flicker 
and motion perception. While only parvocellular lesions had an effect on the monkey’s ability 
to discriminate between subtle color differences, these same lesions did not affect the monkey’s 
ability to detect a single large target whose color differed from its background, even when the 
target was equiluminous with its background. This implies that the M pathway is capable of 
conducting some gross color information and can do so at isoluminance. While lesions in the 
M pathway resulted in a deficit of flicker perception, this was universally true for high temporal 
frequencies only. For low temporal frequencies, lesions in the M pathway had no effect, thus 
demonstrating that the P pathway is capable of transmitting low temporal frequency 
information. (Sekuler and Blake, 1994) 

The P pathway also provides information to bloblike regions of the visual cortex. 
Unlike most ganglion cells, the cells in these bloblike regions are not at all concerned about 
orientation. These cells called blobs exhibit color opponency in a manner similar to that of the 
ON and OFF center cells. However, these cells turn ON and OFF in response to a specific 


chromatic illumination instead of an overall illuminance (Sekuler and Blake, 1994). Having this 


basic physiological background, we can now focus on possible physiological explanations for a 
phenomenon known as the “oblique effect.” 

It is well documented that honzontal and vertical lines are easier to see than lines at 
oblique angles, a phenomenon known as the “oblique effect” (Campbell, Kulikowski and 
Levinson., 1966; Appelle, 1972). This phenomenon has been observed in a variety of visual 
tasks. Essock (1980) divided the oblique effect into two classes. Oblique effects arising from 
basic visual functions such as detection, acuity, and other measures of sensitivity are termed 
Class I oblique effects. The Class I oblique effect is not caused by a bias in the optics of the 
eye (Campbell and Kulikowski, 1966), but is thought to result from the orientation bias of the 
P-cells located in the visual cortex (Lennie, 1974). Several studies have shown that an oblique 
effect is most observable when a stimulus with a high spatial frequency and a low temporal 
frequency is presented to the fovea (Maffer and Fiorentini, 1973; Berkley, Kitterle and 
Watkins, 1975; Camissa, Blake and Lema, 1977). 

Class If oblique effects arise from tasks that require subsequent processing of stimulus 
information. For example, a task measuring the detection threshold of a stimulus onented 
either obliquely or non-obliquely would result in a Class I oblique effect. When an observer 
must press one of two buttons indicating whether two simultaneously-presented stimuli of 
various Orientations are the same, or whether they are different, the result would be a Class I 
oblique effect. This Class Il oblique effect would be the result of encoding or further 
processing of stimulus information required for task completion. An important distinction, 
however, is that the Class I oblique effect discussed above results from achromatic stimuli. For 
chromatic stimuli, the results are not clear and leave uncertainty as to whether an oblique effect 
exists (Kelly, 1975b; Murasagi and Cavanagh, 1989). This thesis explores a Class I oblique 
effect and will be specific as whether this oblique effect is a result of achromatic or chromatic 
stimuli. (Essock, 1980) 

There are various hypotheses explaining why the oblique effect exists. One suggests 
that the world we live in, especially urban areas, contains more stimuli onented honzontally or 
vertically as opposed to obliquely; thus, visual experience plays a role in determining sensitivity 


(Annis and Frost, 1973). However, the oblique effect has been demonstrated with infants as 


1] 


young as six weeks old (Leehy, Moskowitz-Cook, Bnill and Held, 1975). It seems unlikely 
that the visual experience of infants would be sufficient to account for the oblique effect. To 
further confuse the issue, it has been shown that with extensive practice detecting diagonal 
lines, observers may improve their oblique sensitivity until it 1s equal to their sensitivity for 
detecting horizontal and vertical lines (Mayers, 1983; Krebs, 1992). A more widely accepted 
hypothesis suggests that more cells are tuned to vertical and honzontal stimuli than to oblique 
stimuli. 

Similarly disputed is whether the oblique effect is a retinal or a neural phenomenon. 
Campbell and Kulikowski (1966) and Mitchel and Muir (1967) used lasers to bypass the optics 
of the eye by projecting stimuli directly onto the retina. The oblique effect was obtained in 
both studies. Other studies involving head tilt (Rock and Heimer, 1957; Attneave and Reid, 
1968) further investigated whether the oblique effect was a retinal phenomenon. When 
subjects viewed tilted stimuli with their heads tilted the same amount as the stimuli, the stimuli 
were retinally upright, but phenomenally oblique. (“Phenomenally’ descnbes the stimuli’s 
Orientation in the visual frame of reference of the subject [Lasaga and Garner, 1983]. The 
phenomenal frame of reference in this case would be gravitational.) Other, similar studies have 
adopted an arbitrary frame of reference (Rock and Heimer, 1957). In a study by Attneave and 
Reid (1968), subjects were told to think of the top of their heads as vertical, regardless of 
whether or not they were tilted. For head-tilt experiments in which subjects had their heads 
tilted 45 degrees, stimuli that were honzontal/vertical were retinally oblique, but were 
phenomenally upright. Unless otherwise instructed, subjects tended to adopt a phenomenal 
frame of reference rather than a retinal frame of reference. Therefore, when a subject’s head 
was tilted 45 degrees, an oblique effect was obtained for oblique stimuli, even though these 
stimuli were not retinally oblique (Attneave and Olsen, 1967). However, when subjects were 
told to think of the top of their heads as vertical, regardless of their 45 degree head tilt, they 
displayed an oblique effect for stimuli that were gravitationally upright, but retinally oblique 
(Attneave and Reid, 1968). In light of these studies, the oblique effect is highly unlikely to be a 


retinal phenomenon. 


Another possible cause of the oblique effect is a neural phenomenon. This hypothesis, 
and that the origin of this neural phenomenon arises in the P pathway, seem credible even 
without the additional evidence provided by Rabin’s (1992 and 1994) work with Visual 
Evoked Potentials (VEPs). The P pathway is responsible for acuity information and, thus, 
spatial information and is capable of processing low frequency temporal information. The 
Class I oblique effect has been observed primarily at high spatial frequencies and low temporal 
frequencies. 

While the origin of the Class I achromatic oblique effect has been disputed, the dispute 
about a Class I chromatic oblique effect is not over its origin, but over its very existence. One 
of the earliest articles discussing the Class I oblique effect and chromaticity 1s a 1975 study by 
D. H. Kelly. 

The stimuli Kelly used in the experiment were striped luminous-contrast gratings 
flickering sinusoidally. A grating is a pattern of adjacent light and dark bars or stripes, and a 
sinusoidal grating, shown in Figure 2.6, has a gradual transition from light areas to dark areas 


with no sharp edges (Schiffman, 1996). The gratings were presented horizontally, vertically, 
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Figure 2.6. Sinusoidal gratings. Wide stripes correspond to low spatial frequencies. As 
the spatial frequency increases the stripes become thinner. 


and at angles of 45 and 135 degrees (obliquely). A plot of the mean threshold at each 
orientation for various temporal frequencies ranging from approximately 1-40 hertz is shown in 


Figure 2.7. As the figure shows, the thresholds for the oblique presentations are consistently 
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lower than those for the non-oblique presentations for all temporal frequencies. This result is 


consistent with previous results. (Kelly, 1975) 
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Figure 2.7. Achromatic sine-wave flicker sensitivity curves. 10’ = 3 cpd and 22’ = 1.36 
cpd. From Kelly [1975]. 


Although this result was consistent, Kelly wanted to test the hypothesis that luminous 
contrast was a necessary condition for the oblique effect. He repeated the experiment, but 
changed the stimulus to a red-green equiluminous grating. The thresholds obtained, presented 
in Figure 2.8, showed that for the 10-minute stripes, an oblique effect was present for temporal 
frequencies under approximately 10 hertz. However, for the 22-minute stripes, no oblique 
effect was observed. This is due to the luminous grating sensitivity decreasing with decreasing 


spatial frequency, whereas chromatic grating sensitivity remains constant (Kelly, 1975). Kelly 
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Mullen’s article, “The Contrast Sensitivity of Human Colour Vision to Red-Green and 
Blue-Yellow Chromatic Gratings,” described an innovation with her experimental design which 
led to measurements without any type of chromatic aberration. By using a large field size, she 
was able to measure thresholds for low spatial frequencies without the reduction in luminous 
sensitivity shown to occur with spatial frequencies below approximately four cycles/deg. 
Instead of using only one chromatic intensity value for all specific spatial frequencies, as had 
often been done before, Mullen used a number of selected points. This provided more 
accuracy. These factors combined to give the chromatic contrast sensitivity function (CSF) 
obtained for red-green gratings a much different look than previously thought. The CSF still 
had the same basic shape, but the cutoff for high frequencies occurred much earlier at 
approximately 10-12 cycles/deg. (Mullen, 1985) 

The CSF is a method used to describe the visual system’s sensitivity to sinusoidal 
waveforms. Contrast, as defined in a CSF, 1s a relative measure that is computed rather 
than measured. Contrast, the difference between stimuli elements, is formally defined as 
the amplitude of a waveform relative to its mean. Therefore, at a mean luminance level of 
5 cd/m’, a sinusoidal grating with a contrast of 50 percent would have a trough of .25 and 
a peak of .75 cd/m’. This same waveform at a mean luminance level of 500 cd/m* would 
still have a contrast of 50 percent if its peak were at 750 and its trough were at 250 cd/m” 
(Schiffman, 1996). The use of sensitivity (1/threshold contrast) in CSF is similar to 
everyday usage; therefore, a low detection threshold is equivalent to high sensitivity. 
(Schiffman, 1996) 

A CSF for spatial frequency is shown in Figure 2.9. Peak sensitivity is found at 
approximately three cycles per degree (cpd), with approximately 50 cpd being the cutoff 
for high frequency acuity. 

Prior to Mullen’s work, studies with red-green stimuli used frequencies more than 
20 cpd and suggested resolution above 25 cpd (Mullen, 1985). Mullen’s work provided a 
template for further research. It is not by coincidence that Murasagi and Cavanagh’s 1989 
article dealing with the chromatic oblique effect, and using red-green stimuli, chose spatial 


frequencies under the 10-12 cpd cutoff proposed by Mullen. 
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This article by Murasagi and Cavanagh further explored earlier work by Kelly 
regarding the oblique effect for luminous, as well as chromatic, stimuli. Kelly had 


postulated the absence of an oblique effect for chromatic channels if the opponent-color 





pathways for humans, like those in monkeys, were not orientation selective. However, 
research published the same year as Kelly’s article (Poggio, Baker, Mansfield, Sillito and 
Grigg, 1975), as well as additional research a few years later (Michael, 1978), revealed 
that monkeys might possess orientation selectivity in their chromatic channel. The 
possibility of the chromatic channel analyzing orientation independently of the luminous 
channel led the authors to design an experiment to test this possibility in humans by 
determining if an oblique effect obtained with chromatic stimuli differed from that 
obtained with strictly achromatic stimuli. (Murasagi and Cavanagh, 1986) 

To test this possibility, the researchers used a constant temporal frequency of 2 Hz 
and spatial frequencies of 2, 4 and 8 cpd were used. The stimuli were sinusoidal gratings. 
The gratings were presented at oblique (45 and 135 degrees) and non-oblique (0 and 90 
degrees) angles. Axial chromatic aberration was taken into account by having the subjects 
view the stimuli through an achromatizing lens. A revised ascending method of limits was 
used to determine thresholds for both luminance and chromatic stimuli. Since the 
production of an isoluminant stimulus is a non-trivial matter, with isoluminance varying 


slightly from subject to subject, Murasagi and Cavanagh used stimuli in five different areas 


\7 


in the neighborhood of equiluminace. The maximum threshold for the spatiotemporal 
region they were investigating was assumed to occur at equiluminance. They made this 
assumption because a chromatic grating with no luminous component should be the 
hardest to detect, as detection is by color alone rather than by color and luminance. 

A significant main effect of orientation was present for all observers, as was a 
significant three-way interaction between grating types. Spatial frequency and orientation 
effects were present for three of the four observers. This showed that, for three of the 
four observers, the effects of the four orientations at certain spatial frequencies were 
different for achromatic and chromatic stimull. 

Like Kelly, Murasagi and Cavanagh have possible problems with spurious 
luminous components. By taking five measurements in the neighborhood of 
equiluminance, the contribution of these components has probably been reduced. 
However, if the actual isoluminance point were in a region between the areas they picked, 
then a luminous component would be present in their stimuli. Additionally, while spatial 
frequencies of 2, 4 and 8 cpd were used, a chromatic oblique effect was present only at 8 
cpd for three of the four observers. At 8 cpd, chromatic aberration is a factor. The 
researchers used an achromatizing lens to account for axial chromatic aberration, but they 
did not take into account lateral chromatic aberration. Post hoc analysis minimized the 
possibility of this being a factor. As the authors themselves state, any slight misalignment 
between a subject and the achromatizing lens would result in a spurious luminance 
component. (Murasagi and Cavanagh, 1986) 

During the same time frame as the Murasagi and Cavanagh work, Bradley, 
Switkes and De Valois, were also exploring Kelly’s earlier work. The authors designed an 
experiment to compare the visual processing of chromatic and luminance information. 
The prolonged viewing or adaptation of a sinusoidal grating desensitizes the observer to 
similar gratings, especially when the similarity is in onentation or spatial frequency. 
However, this desensitivity is not present for gratings with orientations differing by 
approximately 45 degrees or spatial frequencies differing by 1.5 octaves. Thus, this 


adaptation has been termed selective. 
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The effects of selective adaptation have been used as psychophysical evidence for 
the presence of spatial frequency-selective and orientation-selective neurons in the human 
visual system. However, the behavior of cells displaying selective adaptation for spatial 
frequency when measured psychophysically has not always been consistent with their 
physiology. Thorell, De Valois, and Albrecht (1984) observed neurons that displayed 
different spatial frequency tuning depending on whether the stimulus contained luminance 
or color. They observed low-pass tuning for chromatic gratings, but band-pass tuning for 
luminance gratings. Additional studies by Livingstone and Hubble (1984) and Lennie, 
Sclar and Krauskopf (1985) found that cells in the visual cortex that responded to 
isoluminant color contrast did not display selective adaptation for orientation or spatial 
frequency. (Bradley, Switkes and De Valois, 1988) 

Bradley et al. (1988) set out to explore this inconsistency with spatial frequency 
adaptation and orientation for chromatic gratings. The zero contrast condition for all 
gratings was a uniform yellow field with a chromaticity that was adjusted for each 
observer’s differing sensitivity to red and green phosphor emissions. This varying of the 
zero contrast condition enabled presentation of the red-green sinusoidal gratings at each 
observer’s isoluminance axis. Both isochromatic and isoluminant gratings were presented 
and were viewed through an achromatizing lens. To overcome problems associated with 
making repeated measurements of a decaying effect, the researchers used a long inuitial 
adaptation period followed by alternation of a brief stimulus presentation with a brief 
adaptation period. The stimuli used for adapting was a 2 cpd grating, run separately for 
each of the four possible conditions of horizontal or vertical and luminance or chromatic. 
For a spatial frequency of 2 cpd, thresholds for luminance gratings were similar in both 
pre- and post-adaptation trials for oblique angles and showed the desensitivity expected at 
horizontal and vertical angles. However, while the pre-adaptation data for chromatic 
gratings at this frequency did not show an oblique effect, an oblique effect was evident in 
the post-adaptation data. (Bradley, et al., 1988) 

A similar experiment varying spatial frequency while keeping orientation constant 


confirmed that, for varying spatial frequencies, a specific spatial frequency adaptation 
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effect can be observed for sinusoidal luminous gratings. This experiment was repeated 
with chromatic gratings, and a specific spatial frequency adaptation effect was also 
observed with sinusoidal chromatic gratings. The results of this study along with previous 
psychophysical studies demonstrating the parallels between the data for luminance and 
color (De Valois and Switkes, 1983; Switkes and De Valois, 1983; and Ware and 
Mitchell, 1974) suggest that for the beginning stages of human vision, color and luminance 
are processed in a similar manner. (Bradley, et al., 1988) 

In addition to the numerous psychophysical studies on the oblique effect, other 
studies have been conducted electrophysiologically. When studying primates or other 
animals, collecting data is often not possible through psychophysical means. Although an 
animal may not be able to verbalize or react, a response may still be obtained 
electrophysiologically. By electrically stimulating an individual cell, it is possible to 
monitor the cell output. VEPs provide an additional method of studying the role of 
chromatic patterns in perception (Rabin, 1994). In Rabin’s 1992 paper “VEP’s in Three- 
Dimensional Color Space,” a Class I oblique effect at isoluminance or a chromatic oblique 
effect was shown at the spatial frequency of 1 cpd. Psychophysically, the Class I oblique 
effect for luminance or chromaticity is typically not obtained at low spatial frequencies. 

The Class I oblique effect for achromatic stimuli has been obtained under a number 
of different conditions. It has been demonstrated psychophysically (Campbell and 
Kulikowski, 1966; Camisa, et al., 1977) and electrophysiologically (Maffei and Campbell, 
1970; Rabin, Switkes, Crognale, Schneck and Adams, 1994). Electrophysiologically the 
oblique effect 1s evident by comparing the output of microelectrodes for oblique and non- 
oblique stimuli. These microelectrodes monitor the VEPs of cells as they are exposed to 
different stimuli. Psychophysically, the Class I oblique effect is evident by comparing the 
responses of subjects for oblique and non-oblique stimuli. A Class I chromatic oblique 
effect could be measured the same way. However, the Class I oblique effect for chromatic 
stimuli has not been obtained under various conditions, and whether such an oblique effect 
actually exists is a matter of debate. To participate in this debate, it is necessary to 


understand how the information that the eye collects is processed. One explanation is that 
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the eye uses a process similar to Fourier analysis so named after the nineteenth century 
mathematician responsible for this analysis. 

Jean Baptiste Fourier studied how heat flows through an object when it is heated 
up and found that heat behaved in waves. He modeled these waves using complex 
equations, and discovered that they consist of periodic waveforms. Fourier found that any 
quantity that changed in a complex manner over time could be converted into a series of 
simple sinusoidal functions. Each sinusoidal pattern could be defined by its period, 
frequency, and angular velocity. This process is now known as Fourier analysis. (Who is 
Fourier?, 1995) 

Fourier analysis can be used to analyze a natural scene by decomposing it into a 
sum of a series of sinusoidal components, each having a different spatial frequency, 
amplitude, and orientation. Vision scientists believe that the human visual system uses a 
process similar to Fourier analysis to process visual imagery. The human eye receives 
different intensities of light reflected from an object. These light intensities pass through 
the cornea and filter down into the photoreceptors. The photoreceptors then send an 
electrical signal to the brain, where these neural responses are categorized into specific 
spatial channels. Psychologists believe that this sensory input 1s transformed into a neural 
response, which 1s then categorized into a perceptual experience. If the visual system 
passes the image, and this image corresponds to a perceptual experience, then the observer 
can recognize the object. However, if your cornea is degraded--e.g., a cataract--high 
spatial frequency sinusoidal waves will not pass through the lens and be sent to the brain. 
A degraded signal such as this or a lack of a similar perceptual experience may result in 
failure to recognize the object. The absence of high spatial frequencies will cause the 
image to appear blurry. In some cases, the amplitude of those missing high spatial 
frequencies can be increased so that these signals can be sent to the brain. 

Scene (Figure 2.10) can be broken down into many different visual components by 
Fourier analysis or other tools. These components or parameters include color, 
orientation, and spatial qualities and include the scene as a whole, as well as for the 


individual objects that comprise the scene. By manipulating the parameters of an object 


ZA 


(e.g., armored personnel carrier [APC]) in a scene, it is possible to camouflage this object. 


This may be done by changing the object’s color to match that of its background through 
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Figure 2.10. Armored Personnel Carner 


temporary means such as netting, or by more permanent means such as paint. A Fourier 
analysis shows the spatial composition of the scene. The low spatial frequencies are 
located in the center of Figure 2.11, and the high spatial frequencies are found in the 
comer regions of the figure. The high spatial regions result from the edges of the APC. 
These high spatial frequencies contrast with the low spatial frequencies found elsewhere in 
the scene. Netting would reduce these high spatial frequencies and would also lessen the 
edge effect evident in Figure 2.12, thereby enhancing the APC’s camouflage. 

We have looked at a scene’s color, orientation, and spatial information and how 
these parameters can be manipulated to achieve better camouflage. The parameters 
manipulated in the experiment in this thesis are spatial and orientation. Since the Class I 
oblique effect has been primarily observed at low temporal and high spatial frequencies, a 


temporal frequency of 2 Hz and spatial frequencies of 1,3 and 7 cpd were chosen. The 
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Figure 2.11. Fast Fourier Transformation on the APC. 








Figure 2.12. High-pass filter of the APC originally shown in Figure 2.10. 


red-green spatial CSF begins to decrease at frequencies greater than 1 cpd; therefore the 
frequencies of 3 and 7 were chosen knowingly, trading off sensitivity for the advantages of 
a higher spatial frequency where there would be a higher likelihood of observing an 
oblique effect. Frequencies higher than 7 cpd were not chosen due to increasing effects of 
chromatic aberration. 

A Class I chromatic oblique effect was expected to be observed at spatial 
frequencies of 3 and 7 cpd. Psychophysically, the Class I oblique effect has not been 
readily observed at spatial frequencies as low as 1 cpd and accordingly a Class I chromatic 


oblique effect was not expected to be observed at this spatial frequency. 
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Il. METHODS 


A. SUBJECTS 


The experiment was conducted concurrently at the Naval Postgraduate School 
(NPS) and the University of Louisville, Kentucky (UL). Four subjects, 2 NPS and 2 UL, 
volunteered for this experiment. All subjects had normal (20/20), or corrected to normal, 
acuity and color vision. Color vision was verified with pseudo-isochromatic plates. Two 
of the four subjects (1 NPS and 1UL) were naive as to the purpose of the experiment. 
The other two subjects, the author and the remaining UL subject, were experienced 
psychophysical observers. All subjects signed an informed consent and were briefed on 
the ethical conduct for subject participation specified in the Protection of Human Subjects, 
SECNAV Instruction 3900.39B. Subjects were screened for uncorrected astigmatic 


errors by determining spatial resolution limits for 0°, 45°, 90°, and 135°. 


B. APPARATUS 


Stimuli were presented by a VisionWorks computer graphics system (Vision 
Research Graphics, Inc.) on an IDEK MF-8521 high resolution color monitor (21" X 20" 
of viewable area) equipped with an non-glare, anti-reflect, P-22 phosphor. The monitor 
had a resolution of 800 by 600 pixels (x=75.02 and y=74.92 pixels/degree), 98.9 Hz 
frame-rate, mean chromaticity of r= 0.334 , g = 0.336, b =0.300 (1931 CIE), anda 
maximum luminance of 100 cd/m2. Refer to Table 3.1 for the chromaticity and 
luminance coordinates for each phosphor. The University of Louisville’s apparatus and 
procedure were identical to the Naval Postgraduate School’s, except that the stimuli were 
displayed on a 17” Nanao Flexscan F2.21 color monitor. Subjects viewed the monitor 
from 1.5 meters and were positioned by an adjustable chinrest. A small floor lamp (2.6 


cd/m2) was positioned behind the monitor to reduce screen glare. 
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x y Zz Luminance (cd/m”) 
Red phosphor 617 345 .038 24.0 
Green phosphor .334 .581  .085 88.7 
Blue phosphor 162 081 .757 12.7 


Table 3.1. Chromaticity and luminance of monitor 


STIMULI 


Sinusoidal gratings were presented within a spatially windowed circular test field that 
subtended 7.59° of visual angle. The Gaussian window was truncated at +1 standard 
deviations for both x and y directions. The test patterns were one-dimensional spatio-temporal 
sinusoids of varying orientation (principal and oblique), spatial frequency (1.0, 3.0, and 7.0 
cycles/degree), and color contrast. Test patterns for each subject consisted of two onentations, 
principal (0° and 90°) and oblique (45° and 135°). For each subject, maximum sensitivity for 
each orientation within the principal and oblique grouping was chosen. All sinusoids were 
raised cosines temporally modulated at 2.0 Hz. The sinusoid pattern was presented in a 
1500 msec interval with contrast ramped on and off according to a linear window. 
(Contrast peaked at 202 msec and fell at 1304 msec). 

Color contrast was computed by different ratios of percent red and green (Sellers 
et al.,1986). The monitor was controlled by a Cambridge Research Systems VSG 2/4 
video board that was linearized to 10 bits of resolution per gun. The outputs of each gun 
were linearized by means of stored look-up table file. Sixteen different sinusoidal red- 
green color mixtures were generated by changing the red phosphor only, green phosphor 
only, or by changing the red and green guns in fixed proportions. Color contrast was 
defined according to the (Michelson) formula shown in Equation 3.1. Blue gun was held 
constant in all quandrants. Red and green gun values were used in the determiniation of 


red and green contrast as shown in Equations 3.2 and 3.3. 
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contrast = (peak - trough) / (peak + trough) Bal 
red contrast = (red gun value - 50) / 50 32 
green contrast (green gun value - 50) / 50 33 


D. PROCEDURE 


Thresholds were determined by a two-alternative forced choice adaptive 
psychometric procedure, QUEST (Watson and Pelli, 1983). Threshold was defined at 75 
percent correct. A total of 480 trials, 30 tnals per condition, were randomly presented 
within each session. A session (~ 45 minutes) consisted of one sinusoidal condition with 
16 different red-green color mixtures. “A subject’ had to complete six sessions to 
contribute one threshold point for all conditions. 

At the beginning of each session, subjects dark-adapted for approximately five 
minutes before initiating the first experimental trial. Three of the four subjects were tested 
monocularly, while the fourth subject (UL) was tested binocularly. At the beginning of 
each trial, the subject was instructed to focus on a fixation cross (.19° by .13°) located in 
the center of the screen. The subject initiated the first trial with a keyboard response, the 
fixation cross extinguished followed by presentation of the first interval, 121 msec ISI, and 
then presentation of the second interval. The subject’s task was to detect which interval 
contained the sinusoidal grating. The next trial followed 250 msec after the subject’s 
keyboard response. 

Color contrast thresholds were determined from 16 different color-mixture ratios. 
The sixteen different ratios could be divided into four different percent red and green 
quadrants. Quadrant one started with 100 percent green and 0 percent red, quadrant two 
started at 100 percent red and 0 percent green, quadrant three started at -100 percent 
green and 0 percent red, and quadrant four started at -100 perecent green and 0 percent 
red. The red-green ratios within each quadrant were 0.0, 0.5, 1.0, and 2.0. The 


thresholds from these red-green ratios will form an ellipse with the half-length of the axis 


vay 


measuring color discrimination and the half-width of the axis measuring brightness 


discrimination (Sellers, 1986). 
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IV. RESULTS 


Many researchers have carried out experiments to determine visual sensitivity to color 
differences. One way of determining these differing sensitivities is through differing values of 
International Commission on Illuminance (ICI or, more commonly, CIE, for the French translation) 
primaries (Kaiser and Boynton, 1996). These primaries allow all colors to be specified in terms of 
three numbers representing the red, green and blue primaries. The color-matching type experiment 
is set up to test if an observer can discriminate between a chosen color and another color similar to 
this color. Color-matching experiments have looked at the standard deviations of color-matchings 
for representative colors throughout the color spectrum. These standard deviations are directly 
related to the corresponding just-noticeable difference of colors. (Brown and MacAdam, 1949) 

Numerous surveys of differential thresholds have been carried out, but W. D.Wright (1941) 
and D.L. MacAdam (1942) completed two of the more extensive surveys. “MacAdam plotted the 
results of the survey on the chromaticity diagram in terms of the standard deviation of color-matching 
in several directions for selected colors.”(Brown and MacAdam, 1949) The figures resulting from 
this survey formed closed curves on the diagram and the closed curves were elliptical in shape. 
(Brown and MacAdam, 1949) 

In a later paper, Silberstein and MacAdam discussed that the errors of these closed curves 
were Normally distributed. They deduced that the curves should be ellipses, as they appeared to be. 
They further expostulated that if the variations were not confined to chromaticity, the closed curves 
would form ellipsoids rather than ellipses. Using the assumption that the probability of making a 
match that falls within a specific region of color space, near the target color, was not changed by any 
change of primaries, Silberstein proved the standard deviation figures to be ellipsoids. (Brown and 
MacAdam, 1949) 

The fact that there was a theoretical explanation for MacAdam’s ellipses and not just an 
empirical observation was interesting. However, the discrimination ellipsoids of Brown and 
MacAdam were obtained for a bipartite field only. Noorlander, Heuts and Koenderink (1979 and 
1980) and Noorlander and Koenderink (1983) furthered the work of Brown and MacAdam by 


ue 


extending their methodology from a simple bipartite field to more complicated stimuli by varying 
temporal and spatial frequencies. When three primaries are used discrimination ellipsoids are 
obtained. However, if just two primaries are used, such as red and green, than a cross-section of a 
discrimination ellipse is obtained. Neurlander, Heuts and Koenderink (1980) did this by obtaining 
discrimination ellipses for a number of different spatial and temporal frequencies. The lengths of the 
major and minor axes of these ellipses, as well as their orientation, are highly dependent on both the 
spatial and temporal frequencies. 

The finding that the cross-section of a discrimination ellipsoid is an ellipse was also used by 
Sellers et al. (1986) in a study of congenital and acquired color defects. However, in the Sellers et 
al. paper, the axes of their graphs were not primaries as in Brown and MacAdam’s or Neurlander , 
Heuts and Koenderink’s, but were percent red’contrast and percent green contrast. Even without 
knowledge of discrimination ellipsoids, the fact that Sellers et al.’s data form an ellipse can be 
explained by examining the following model. The central dashed line in Figure 4.1 is the 
equiluminance axis. Assume there are two luminance processes for detecting the brightening and 
darkening of a spot. The thresholds of these processes are displayed in Figure 4.1 with dashed lines 
and are labeled “BRIGHT” and “DARK.” Similarly, the processes for detecting color are labeled 
“RED1” and “GREEN1.” (“RED2" and “GREEN2" refer to two different thresholds.) As a first 
approximation, the visual threshold will be determined by whichever process has the lowest threshold. 
Therefore, this will be the parallelogram bounded by the four lines “BRIGHT,” “DARK,” “RED1," 
AND “GREEN1.” A phenomenon known as probability summation accounts for the rounding of the 
comers of the parallelogram, and an ellipse is formed. Probability summation occurs near the comers 
of the parallelogram and can be thought of as a sum of two processes, e.g. “BRIGHT” and “RED1." 
For example, if the probability of either of the processes detecting a stimulus is .5, then if both 
processes are independent, the probability of either of them (or both) detecting a stimuli is .75. Thus, 
a contour connecting points where the probability of detection is .5 will exclude the comers of the 
parallelogram (Graham, 1989). (Sellers, 1986) 

A ratio can be determined by dividing the length of the major axis by the length of the minor 
axis. For ellipses with a ratio greater than four, the major axis nearly coincides with the 


equiluminance axis. The angle of discrepancy between the equiluminous axis and the major axis 
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Figure 4.1. Detection model. From Sellers et al.[1986]. 
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The subjects run in this expenment did not have any color defects. Thus, collecting data on 
these subjects for classification purposes was not an exceptionally interesting endeavor. However, 
the possibility of using this methodology to explore a chromatic oblique effect was interesting. If 
oblique and non-oblique sensitivities are the same, and if oblique and non-oblique information is 
processed in an identical manner, then the ellipse obtained from a subject responding to non-oblique 
(horizontal or vertical) chromatic stimuli and the ellipse obtained from the same subject at the same 
temporal and spatial parameters, but with oblique chromatic stimuli, should be identical. A 
“spurious” luminous component is not a problem. Since the data points theoretically form an ellipse, 
the requirement for a point exactly at isoluminance no longer holds. 

The elliptical nature of the data has been used to fit ellipses to the data by the method of 
maximum likelihood. A well-known result from linear regression informs us that the method of 
maximum likelihood 1s identical to the method of least squares in this case (Larsen and Marx, 1986). 
The programs used here to fit ellipses minimize the sum of the squared error. Assistant Professor 
Professor Samuel Buttrey of the Naval Postgraduate School in Monterey, CA. created these 
programs, their sub-programs and other programs of use. He created these programs, which are 


found in Appendix A, using the statistical package S-Plus. 


The following terminology will be used to describe ellipses and their parameters (Figure 4.2). 


terminology definition 
a half-length of the major axis of an ellipse 
b half-length of the minor axis of an ellipse 
0 angle as measured from the x-axis to the major axis of an ellipse 
x, x-coordinate for data point 1. X-axis is percent red contrast 
y, y-coordinate for data point 1. Y-axis is percent green contrast 
ee x-coordinate for the center of an ellipse 
y, y-coordinate for the center of an ellipse 
y distance from (x,, y,.) to a point on the ellipse along the ray 
E, error term 
t polar angle of 7, 
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Five-Parameter Ellipse 





Figure 4.2. Five parameter ellipse. Created by Professor Samuel Buttrey. 


The parameters for the true ellipse are unknown, but they may be estimated from the data. Three 
models were used to represent the underlying ellipse, and ellipses will be classified according to the 
method in which they were modeled. The polar angle relative to the coordinate ¢, is fixed at the start 
of the experiment and is determined by tan(z) = y,x,. The sixteen values for 7, in degrees are (0, +30, 
+45, +60, +90, +120, +135, +150,180). Each of theses ellipses possesses an equiluminous axis along 
which, by definition, luminance is constant. The exact determination of this equiluminous axis is 
difficult, but for ellipses with an a/b ratio greater than four, this equiluminous axis is closely 
approximated by the major axis of the ellipse (Sellers, et al., 1986). A ray that coincides with the 
equiluminous axis will vary in color along 7,, but will have a constant luminosity. Any ray that 1s not 


aligned with the equiluminous axis will have a constant red to green percent contrast ratio along the 


oS 


ray, but luminosity will not be constant along the ray. Both x and y values can be positive or 


negative. 
A. CLASS I-TYPE ELLIPSES 


The Class I-type ellipse has five estimated parameters. These parameters are a, b, 0, x, and 

y,. The predicted r ,or 7,, is a function of the estimated parameters and F, = f(a,b,0, x .y .t). The 
actual model used is r,;= 7, + &. Here & are independent identically distributed (iid) and N(0,o’). 
The function f{a,b,0, x, y, ¢) is complicated and can be found in Appendix A (ell.pred). The 
objective function is the sum of the squared differences between the observed 7, and the predicted 
?.. Data were collected from both oblique and non-oblique stimuli, the class I-type ellipse is an 
ellipse that is obtained by fitting an ellipse to all of the data for a specific spatial frequency. For 
example, Subject One completed five runs at each condition. Each run results in the calculation of 
a threshold along each ray; thus, 16 thresholds were determined for this subject on five different 
sessions. This resulted in 80 data points for the non-oblique condition and 80 data points for the 
oblique condition for each of the three spatial frequencies used. Class I-type ellipses are fitted to the 
combined data of a subject at a specific spatial frequency. For Subject One, 160 data points were 
used, and an ellipse was fitted by the method of maximum likelihood. In the past, the ellipses 
obtained in this manner have been forced to have their center at the ongin (Sellers, et al., 1986). 
However, much better fitting ellipses are obtained by allowing the center not to be pinned to the 
origin. The centers obtained for most subjects were generally close to the origin. The fact that a 
better fit was obtained by letting the ellipse be centered at coordinates other than the origin may be 
an indicator that centers may have some sort of bivariate distribution across subjects, or it may be an 
effect caused by the monitor. However, data from UL were definitely not centered on the origin and 
tended to have a center in the second quadrant. The ellipse programs in Appendix A allow the ellipse 
center to be pinned at the origin (fit.center=F) or to “float” (fit.center=T), but it did not make sense 
to pin the center to the origin when some of the actual ellipse centers obtained were definitely not at 


the origin. 
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B. CLASS-I TYPE AND CLASS-Iil TYPE ELLIPSES 


The Class-I Type ellipse has six (program ellipse.II) and the Class-III Type ellipse has seven 
(program ellipse.III) estimated parameters. The models for the 7, for the three classes of ellipses are 
shown below. The model for Class-I Type ellipses may be changed so that it is ffa, b + 6,, 9, x, 
and y.) by changing which.type from one to two in the program ellipse.II. 


F (class I type ellipses) = fla,b,0,x.,y.) 4] 
F(class II type ellipses) = fla+6,,b,8,x,,y.) 42 
F(class III type ellipses) = fla+6,,b+6,,0,x,,y.) 43 


The Class-II Type and Class-III Type ellipses use the information of whether the data were 
from an oblique or non-oblique condition. An ellipse is then fitted to the data, but the additional 
information of whether the data were from an oblique condition or a non-oblique condition is used 
to determine a 6, and/or a 6,. This is done by fitting a Class-II or Class-III Type ellipse to the data. 
Actually, two ellipses are fit to the data, one to the non-oblique data and one to the oblique data; but 
a common ellipse center and theta are maintained for both ellipses. If the sum of the objective 
functions for the ellipse fitted to the non-oblique data and the ellipse fitted to the oblique data is 
smaller than the objective function for the ellipse fitted to all of the data, then there will be a 6, and/or 
aO,. Ifthe 6, or 6, are small (This will be quantified in the next section.), then this orientation 
information does not significantly improve the fit of the ellipse. Ifthe 6, and 6, are large, then this 


additional information significantly improves the fit of the ellipse. 


C. STATISTICAL TESTS 


A complete discussion of the statistical test used for comparing ellipses can be found on pages 
103-104 of Bates and Watts (1988). This test is an approximation, due to non-linearity, of an F-test. 
The derivation of this F statistic and how it is obtained are shown in Table 4.1 and Equation 4.4 


respectively. For this experiment, the full model consists of Class-III type ellipses, and the partial 
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model consists of Class-II type ellipses. The full and partial models were chosen this way because 
of interest in the significance of the difference in the length of the major axis between non-oblique and 


oblique data. 


Sum of Degrees of Mean Square _ F Ratio 
Source Squares Freedom 


Extra parameters S,=S,-S, v,=P;- P, s =S./v, SE Se 
Full model v-=N-P, s° = S,/ v, 


Partial Model 
Table 4.1. Extra sum of squares analysis. From Bates and Watts [1988]. 





S/v 
é”’e approximately a an 
Sv, ~ aaa 


If the equiluminous axis 1s identical to the major axis of an ellipse, then if the length of the 
major axis for oblique stimuli is greater than the length of the major axis for non-oblique stimuli, a 
chromatic oblique effect has been observed. Additionally, if the length of the ellipse minor axis for 
oblique stimuli is greater than the length of the minor axis for non-oblique stimuli, a luminous oblique 
effect has been observed. For ellipses with a major axis to minor axis ratio of four or greater, the 
equiluminous axis is closely approximated by the major axis (Sellers et al., 1986). The ellipses 
obtained generally had a ratio less than four, but the major axis was still used as the equiluminous axis 
as a rough approximation. 

Four subjects took part in this experiment. Data from Subjects One and Three were collected 
at the NPS, while data from Subjects Two and Four were collected at UL. The data for these 
subjects at a spatial frequency of one cpd and both oblique and non-oblique orientations are shown 
in Figure 4.3. The measurements for a, b, theta, x, and y, are shown on the individual graphs and are 


displayed in Table 4.2, as well. 
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. . Sub 1 1 cpd 2 hz . Sub 2 1cpd 2 hz 
non-oblique = points/solid line —_ oblique = pluses/dashed line non-oblique = points/solid line oblique = pluses/dashed line 


percent green contest 
percent green contratt 





percent red contrast percent red contrast 


non-oblique e® oe? b= 0.03424 oblique as 0.04691 ila non-oblique a= 0.05541 be 0.02196 oblique a= 0.04845 be 0.01964 
thete= -1.527 x« 0.0045 y= 0.0006 p-value= 0.975 theta= -11.349 x= -0.0028 y= 0.0057 p-value= 0.091 
Sub 3 1 cpd 2 hz Sub 4 1 cpd 2 hz 


non-oblique = points/solid line —_ oblique = pluses/dashed line non-oblique = points/solid line oblique = pluses/dashed line 


percent green contact 
percent green contast 





-0.10 -0.05 0.0 0.05 0.10 -0.10 -0.05 0.0 0.05 0.10 
percent red contrast 
; : percent red contrast 
Ty ae geht! Again a non-ablique = 0.09675 b= 0.02234 cbique e= 0.09081 b= 0.0306 
theta= 36.733 x= 0.0059 y= -0. 2 theta= -15.689 x= -0.0033 y= 0.0053 p-valuee 0.513 


Figure 4.3. Ellipses and values for Subjects 1-4 at 1 cpd 









Ellipse Type |Subj|__a__ | __b__|_ theta_|center.x|center.y| delta.a | delta.b | *p-value| 
nomblgque | 1 {04703 asagt|~1.5274 [0.004«7 0.00064 [| 
eblique | 1_| 0.04691 | 0.0351 | -1.5274 | 0.00447 | 0.00064 | -0.0001 | 0.0009 | 0.975 | 
non-oblique | 2 | 0.05541 | 0.02196 |-11.349 | -0.0028 |o.005e9, | S| 
non-oblique | 3 | 0.06444 | 0.04976 | 36.733 |0.00593/ -0.004 | | | 
Inon-oblique | 4 | 0.09675 | 0.02294 | -15.689 | -0.0033/0.00525; | | 
he EE es ee ee ee ee 
ee eee ee | ee eee 


Table 4.2. Ellipse values for Subjects 1-4 at 1 cpd 












oe 


The p-values were obtained by the F-test approximation discussed in the models section; 
however, an example is shown below for further clarification. For Subject One the ellipse from the 
non-oblique data and the ellipse from the oblique data are almost identical, whereas for Subjects Two 
and Three the opposite of a chromatic oblique effect (oblique ellipses with shorter a’s than non- 
oblique ellipses) is displayed and for Subject Four, a chromatic oblique effect is shown. However, 
with an alpha of .05, none of these results are significant. 

Here is an example of how the p-values were calculated. Subject 1 completed 5 runs at all 
conditions. For each run, a total of 16 thresholds were calculated, so for a spatial frequency of 1 cpd 
a total of 80 oblique data points and 80 non-oblique data points were collected for this subject. The 
data was input to the program ellipse. IIL, and the ellipse center was allowed to float or not be pinned 
to the origin. From the output of this program, the objective function value is obtained. This 
objective function value is the sum of the squared error (SS,, for Class-IIJ Type ellipses). This value 
is subtracted from the objective function value obtained from the output of a Class-II Type ellipse 
obtained with the same data, SS,. The difference is then divided by the difference in the number of 
estimated parameters, or degrees of freedom, between the two classes of ellipses. This is the 
numerator for the equation. There is only one additional estimated parameter for Class-III Type 
ellipses, compared to Class-II Type ellipses, so this number is a one. Finally, the denominator is the 
value of the objective function obtained from the Class-III Type ellipse output divided by its degrees 
of freedom. For Class-JJI Type ellipses, seven parameters are estimated, so the 160 degrees of 
freedom for subject One decreased to 153 degrees of freedom. The resulting fraction 1s shown in 
equation 4.5. This fraction is referred to the F distribution with 1 and 153 degrees of freedom. Thus 


the fraction is approximately an /' random vaniable with degrees of freedom given by 1 and 153, 


(SS,, - SS)/Q) 
SS,,,/153 


4.5 


if the hypothesis that the ellipses differ only by 6,, and if iid Normal errors are true. A p-value is then 


calculated through tables or statistical programs. P-values for other subjects were calculated 
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similarly, with the degrees of freedom reflecting the number of observations the subject had for that 
condition. 

The data for subjects at a spatial frequency of three cpd and both oblique and orientations are 
shown in Figure 4.4. The measurements for a, b, theta, x, and y, are shown on the individual graphs 
and are displayed in Table 4.3, as well. A chromatic oblique effect is shown for Subjects One, Two 
and Three and is significant for Subjects One and Two. 

The data for the subjects at a spatial frequency of seven cpd and both oblique and non-oblique 
orientations are shown in Figure 4.5. The measurements for a, b, theta, x, and y, are shown on the 
individual graphs, and are displayed in Table 4.4. An achromatic oblique effect is expected here and 
is evidenced by 6,<>0 leading to larger 5 values for oblique ellipses compared to the b values for non- 
oblique ellipses. A peculiarity of the program that determines the 6 values is that if the 6 value is 
negative and if it is larger in magnitude than the value to be added to, then the signs of both values 
must be reversed. The hypothesis of 6,<>0 was tested in a manner similar to 6, <>0, and all subjects 
displayed an achromatic oblique effect. This achromatic oblique effect was significant for Subjects 
One, Two and Four. A chromatic oblique effect is shown for Subjects One and Two, but is not 
significant for either. 

The data collected from the subjects were extremely variable. This variability is not only from 
subject to subject, but also from day to day and run to run. To display some of this variability, Table 
4.5 shows p-values (uncorrected for multiple comparisons) for a subject’s run at a specific condition 
against all of the other runs at this identical condition. 

In summary, at one cpd neither an achromatic nor a chromatic oblique effect was shown. At 
three cpd a chromatic oblique effect was shown for three subjects and was significant for two of 
them. At seven cpd both achromatic and chromatic oblique effects were shown. All four subjects 
showed an achromatic oblique effect and this oblique effect was significant for three of them. Only 
two of the four subjects showed a chromatic oblique effect, but neither of the p-values were 


significant. 
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Sub 1 3cpd 2 hz Sub 2 3cpd 2 hz 
non-oblique = points/solid line —_ oblique = pluses/dashed line non-oblique = points/solid line — oblique = pluses/dashed line 


percent gteen contvast 
percent green contast 





-0.10 -0.05 0.0 0.05 0.10 -0.10 -0.05 0.0 0.05 0.10 
percent red contrast percent red contrast 
non-oblique = 0.05594 b= 0.02284 oblique e* 0.06445 b= 0.02256 non-oblique a 0.0459 b= 0.01365 oblique e= 0.0556 b= 0.01554 
theta= - 19.535 x= 0.0122 y= 0.0001 p-value= 0.028 theta= -16.952 x= -0.0153 y= 0.0131 p-value= 0.005 
Sub3 3cpd 2hz Sub 4 3cpd 2 hz 
non-oblique = points/solid line oblique = pluses/dashed line non-oblique = points/solid line _ oblique = pluses/dashed line 


percent green contest 
percent green contrast 





-0.10 -0.05 00 0.05 0.10 0.2 -0.1 0.0 0.1 02 
percent red contrast 
non-odlique e= 0.1722 b= 0.02344 oblique e= 0.20257 b= 0.03495 
thets= -17.602 x= -0.0054 y= 0.0093 p-value= 0.268 


percent red contrast 
non-oblique a= 0.08941 b= 0.02025 oblique e= 0.10976 b= 0.02385 
theta= -12315 x= 0.008 y= 0.0023 p-value= 0.057 


Figure 4.4. Ellipses and values for Subjects 1-4 at 3 cpd 


Ellipse Type |Subj| a__ | __b__—|_—ttheta_| center.x|center.y| delta.a | delta.b | *p-value| 
non-obique [1 W0S554 Ouzxe4 Oszs5 ootaee|cooot| PD 
oblique | 1_—| 0.06445 | 0.02256 | -0.3235 | 0.01224 | 0.0001 | 0.0085 | -0.0003| 0.028 | 
non-oblique | 2 | 0.0459 | 0.01365 | -0.2959 | -0.0153 |0.01315| | | 
oblique | 2 | (0.0556 | 0.01554 | -0.2959 | -0.0153 | 0.01315 | 0.0097 | 0.0019 | 0.005 | 
non-oblique __3_| 0.08941 | 0.02025 | -0.2149 | 0.00796 |0.00232/ | | 


non-oblique 4 | 0.1722 | 0.02344 1-0.3072|-0.0054| 0.0093, | | 
FHodeltaa=0] 1] a 
[Ha delta.ao0) | SSCS 





Table 4.3. Ellipses and values for Subjects 1-4 at 3 cpd 
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percent green contrast 


Sub 1 7cpd 2hz . _ Sub2 7epd 2 hz 
non-oblique = points/solid line oblique = pluses/dashed line non-oblique = points/solid line —_ oblique = pluses/dashed line 


A 


percent green contrast 





0.10 0.05 0.0 0.05 0.40 


percent red contrast 
non-cblique e# 0.07852 b= 0.01751 ocblique a= 0.09241 be 0.02723 


percent red contrast 
non-obliqie 2* 0.06371 b= 0.02391 oblique e= 0.06646 b= 0.03584 


theta = - 20.009 x= 0.0052 y= 0.0001 p-value= 0.588 theta= -16,684 x= -0.006 y= 0.008 p-value= 0.067 
Sub 3 7cpd 2 hz | _ Sub4 7epd 2 hz 
non-oblique = points/solid line oblique = pluses/dashed line non-oblique = points/solid line —_ oblique = pluses/dashed line 





percent green contrast 





percent grean cantrast 





-0.10 -0.05 0.0 0.05 0.10 °0.2 “0.1 0.0 0.1 0.2 


percent red contrast 
non-oblique 2 0.21246 b=0.04223 abhque e=0.1789 b= 0.04969 
thetas -13.456 x= 0.0081 y= 0.0018 p-value= 0.43 


percent red contrast 
non-oblique e= 5.35208 b= 0.0178 obliqe e= 0.11977 b= 0.02158 
theta= -14.176 x= -0,0077 y= 0.0024 p-value= 0.059 


Figure 4.5. Ellipses and values for Subjects 1-4 at 7 cpd 





Ellipse Type [Subj] a |b | theta |centerx|centery| delta.a | delta |"p-value 
non-oblique | 1 | 006371] 0.02501) 03403] 0.00817, 0.00012] 
oblique | 1 | 0.06646) 0.03584| -0.3492| 0.00517) 0.00012| 0.0028] 0.0119] 0.588 | 
non-oblique [2 | 0.07862[ 0.01751| 0.2812 -0.006| 0.00788 ff 

oblique | 2 | 0.09241] 0.02723] -0.2912| -0.006] 0.00799| 0.0139| -0.0447| 0.087 
non-oblique | 3 | 6.35208] 0.0178] 0.2474) 0.0077, 0.00238, 
oblique | _3_| 0.11978] 0.02158| -0.2474) -0.0077| 0.00238| -5.2323| -0.0304) 0.059 _ 
non-oblique | 4 | 0.21246) 0.04223) -0.2349| 0.00806/0.00179) | | 
time == [| =) leet oa 

tree ee 


Table 4.4. Ellipses and values for Subjects 1-4 at 7 cpd 













4] 










dL SSCS 
[run [Ron-oblique]_ oblique [non-oblique] oblique |non-oblique] oblique _ 
rf 2 | 0.056} 0.199] 0.242{ 0.048] 0.191] 0.393 
3 __ oat} 0.544) 0.680{ 0.000] 0.558) 0.933 
[4 __ 0399} 0.673] 0.152 0.181] 0.529] 0.001 
Tf s_f__o.066efo.ast_0.722{ 0.141] 0.762] 0.9751 
ee, 9) ee ee) ee 
[| 2 | 0.070} 0.800 0.005] 0.620] 0.899] 0.146 
Cf 3 | 0.222) 0.019 0.900{ 0.014) 0.462) 0.013 
[4 [-__ 0.002, 0.004] 0.002] 0.676 0.000] 0.126 
ee | oe) eo ) eee | oe 
r {| 2 | o770of 0.033] 0.796) 0.095] 0.006 0.231 
[| 3 | 0226 0.017] 0.079] 0.001] 0.101] 0.004 
len) 6d)! ee | 
| 2 | _o7r7| 0.258] 0.546] 0.042 _0.513[ 0.000 
| 3 | 0397 0.926] 0.415] 0.000 0.136] _0.000 
ft 4 dT 620 va] sna 0.659 0.003 nda 





Table 4.5. P-values for comparing one subjects run against their other runs at that same condition. 
P-values have not been corrected for multiple comparisons. Bold values are less than .05/(number of 
runs at that condition). 
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V. CONCLUSIONS 


The purpose of this thesis was to investigate the human visual system’s ability to 
detect certain simple targets. This thesis investigated how an object’s spatial, temporal, and 
color features affected humans’ detection of objects. The results showed that certain spatial 
and chromatic qualities do, indeed, inhibit detection. While real-world objects are much more 
complex than laboratory stimuli, knowledge of spatial and chromatic qualities that inhibit 
detection will assist military designers in the quest for better camouflage. 

Other studies of use to military designers include the numerous studies documenting 
an achromatic Class-I oblique effect and the fact that it is generally found psychophysically 
only at high spatial frequencies. This study produced similar results with all subjects 
displaying an achromatic Class-I oblique effect (p-values of 0, 0, .139 and 0) at a spatial 
frequency of 7 cpd. Previous studies documenting a chromatic Class-I oblique effect, or lack 
thereof, are less useful due to conflicting results and possible problems with luminance 
artifacts tainting results (Kelly, 1976; Murasagi and Cavanagh, 1988). Indeed, the work done 
by Kelly and Murasagi and Cavanagh highlighted the problems in determining a chromatic 
oblique effect due to the difficulty of obtaining isoluminance for a subject. Any deviation can 
lead to the introduction of luminance artifacts and can corrupt the results of the experiment. 
The methodology used in this thesis takes advantage of the elliptical shape of the curve 
connecting thresholds at a fixed temporal and spatial frequency, and makes the exact 
determination of isoluminance unnecessary. 

This thesis supports the hypothesis that a Class-I chromatic oblique does exist. At a 
spatial frequency of 3 cpd, a chromatic oblique effect is evident. A chromatic oblique effect 
is shown for three of the four subjects and, with an alpha of .05, is significant for two of the 
four subjects. Additionally, while the p-value for Subject Three (.057) is not less than .05, 
this subject conducted only three runs with 12 thresholds. An additional run would likely 


reduce this p-value to a value less than .05. 
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The main value of this study is the tool it provides for further investigation of a Class-I 
chromatic oblique effect without the problems associated with luminance artifacts. However, 
this tool does have its drawbacks. The data collection required is extremely time-intensive. 
A total of at least five runs at each spatial frequency is desirable due to the variability of the 
data. With each session lasting approximately 45 minutes to an hour and one run consisting 
of an oblique session and a non-oblique session, the time required for five runs is 
approximately 7.5 to ten hours. This time does not account for three sessions needed to 
determine a subject’s maximum oblique and non-oblique sensitivity. This large time 
commitment on the part of subjects poses problems, as their motivation begins to wane. 
Motivation was a possible problem, as UL subjects did not run consistently; their average run 
time (days) was 19 and 18 versus 10 and 12 for the NPS subjects. This undoubtedly affected 
the variability of their data, as evidenced by the number of their p-values in Table 4.5 less than 
03: 

While the determination of exact isoluminance is not required, it is desirable to 
determine an approximate isoluminance axis. If an ellipse major to minor axis ratio is four 
or greater, then Sellers et al. (1986) state that the major axis is a good approximation of the 
equiluminous axis. Major to minor axis ratios in this study averaged 3.3 and exceeded four 
only one-third of the time. For Subject One at 7 cpd graphically (Figure 4.3), it appears that 
a Class-I chromatic oblique effect is evident. Taking into account that the highest major to 
minor axis ratio for this subject at this spatial frequency is 2.66, a chromatic oblique effect is 
even more likely since the equiluminous axis 1s not likely approximated that well by the major 
axis. In this case, the p-value of .588 does not provide much information regarding what 
actually occurred. The statistical test resulting from the model formulation tests the 
significance in the difference of the length of the axes and does not account for the fact that 
the true isoluminance axis may not be aligned with the major axis. 

This thesis has provided an excellent tool for further research. Possible improvements 
include determination of a subject’s equiluminous axis prior to running the experiment, thus 
enabling the experimenter to choose ratios that would run through this approximate axis. A 
spatial frequency of 3 cpd is worthy of further study with more runs and fewer confounding 


variables (different monitors). 
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APPENDIX A. S-PLUS CODE 


> ellipse 
function(x, y, a = 2, b = 3, e = 0, theta = 0, fit.center = F, grad =T, 
ise enere.hess = T, plot.it = T, chat = T) 


{ 
# This is for a Class-I type ellipse 
# Fit a least-squares ellipse centered at 0 with semi-axes (a, b) 
# and angle to the origin theta, to the data in x, y. The ellipse 
# is here parameterized by a, e (the eccentricity) and theta, 
# in that order. a is always > b. 
# a,b,e and theta are starting points 
# fit.center=F pins the ellipse to the origin when this is true the 
# ellipse is allowed to “float” 
# grad=gradient 
# .hess=hessian 
# plot.it activates plot 
# chat=T shows values as they are computed 
# 
# If ais supplied, and it's a vector, then we've been given 
# starting points for all the parameters. Use ‘em, first making 
# sure that there is the right number (3 if we're not fitting 
# the center, and 5 if we are). 
# 
if(!missing(a) && length(a) > 1) { 
1f((fit.center && length(a) '= 5) || (!fit.center && 
length(a) != 3)) 


stop (paste("Parameter vector has length", length(a), 
", expecting ", ifelse(fit.center, 5, 3))) 


if(length(names(a[2])) == || names(a[2]}) == "e™) { 
if(length(names(a[2])) == 0) 
warning("No param names: using e in pos. 2") 


e <- a[2] 
b <- a[(l1] * sqrt(l - e*2) 


else { 
b <- a[2] 
€ <- sgqrtt(l - (b/a[l])%*2) 
} 
theta <- a[3] 
Li (fl tecenter)) 4 
center.x <- a[4] 
center.y <- a[5] 
} 
a <- a[l] 
} 
else { 
if (missing(a)) 
a) aa oe wader (range (x):) 
if (missing (b) ) 
baa Olmos tt (range (y) ) 
e <- ifelse(a > b, sqrt(l - (b/a)*2), sqrt(l - (a/b)%2)) 
thmimesoing(Eheta))  { 
is2cuee.- ISELtE (x, y) 
theta <- atan(ls.outScoef[2}) 
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} 


center.x <- mean (x) 
center.y <- mean(y) 
} 
Geac= ell Cetx, Vy) 
ibieGe epee okay) 
Graph e<— dev list () (van. oraph ii 
if(length(graph) == || all(is.-na(graph) ) ) 
win.graph () 
else dev.set (graph) 
} 
if (grad) 
grad-funca<— ell orad 
else grad.func <- NULL 
me (fi t.cente rim 


Sstart.vec <- c({a =a, e =e, theta = theta, center.x = 
GenGen.. comer e rs s-sCente aa s 
Hower. Vec <— c(G2ecel, OC-000l, —-2°*" pi, — Int ae 


Upper. Veco see (ini, Cnee2oserm ze pi, Inf, Ing) 
} 
else { 
start.vec <- c{a =a, e =e, theta = theta) 
Lowe havece<—.c(U.0C0l, TOS0COl, > —2e* 1) 
Upper viecew.4— 1G (Uni, 0.999 970, 29% p17) 
} 
out <- nlminb(start = start.vec, objective = ell.res, gradient = 
grad.func, hessian = is.there.hess, lower = lower.vec, upper 
= Upper.-Vec,ate. = tt, mys = <x) myo = Vo loc. te 
plot.it, chat = chat, is.there.hess = is.there.hess, 
fit.centér = fit.center, step.min = 100 * 
-MachineSdouble.eps, scale.upd = 1) # 
p.names <- names (outSparameters) 
cat("In ellipse, check p.names\n") 
browser () 
b <- outSparameters["a"] * sqrt(l - outSparameters["e"]%2) 
if(length (outSparameters) > 3) { 
outSparameters <- c(outS$parameters["a"], b = b, out$ 
parameters [3:length(out$parameters) }) # 
names (outSparameters) <- c("a", "b", "theta", "center.x", 
"center.y") # Beats me... 
} 
else { 
outS$parameters <- c(outSparameters["a"], b = b, outs 
parameters [3}) # 
names (outSparameters) <- c("a", "b", "theta"™) 
} 
outSsigma.sq <- outSobj/length (x) 
outSsigma <- sqrt (outSsigma.sq) 
outSaux <- NULL 
OuESX <= =x 
OuUtSVe<— vy 
if(length(outShessian) > 0) { 
if(qr(outShessian)$rank < ncol(outShessian) ) 
outScov <- "Can't invert Hessian" 
else outScov <- outSsigma.sgq * solve(outShessian) 
} 
OUutLStE <— te 
out$fitted.tt <- ell.tt(x - center.x, y - center.y) 
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pred <- ell.pred(out$Sfitted.tt, outSparameters["a"], 
outSparameters["b"], outSparameters["theta"], 
return.unrotated.too = F,fit.center = fit.center, 
center.x = ifelse(fit.center, out$ 
parameters[“center.x"], 0), center.y = ifelse(fit.center, 
outS$ parameters["center.y"], 0)) 
outSfitted.x <- predSx 
QUES EL ELCdsy = — prErcacy 
outsSftitted.r <= sqrt(predsSx*2 + predsSy%*2Z) 
class(out) <- "ellipse" 
return (out) 


> ellipse.II 

function(x, y, a = 2, b = 3, theta = 0, delta = 0, fit.center 
= TT, 1is.there.hess = T, plot.it = T, chat = T, class.I 
length(x)), which.type = 1) 


FE, grad 
rep (T, 


{ 
# This is for a Class-II type ellipse 
# class.I seperates oblique and non-oblique data. For an x vector of 
# length 160 where the first 80 points were non-oblique data the 
# class.I vector should consist of a boolean vector or length 160 # 

comprised of 80 T’s followed by 80 F’s 
# which.type=1 when testing differences in the major axes (a’s or # 

chromaticity) 
# which.type—-2Z when testing differences in the minor axes (b’s 
# or luminance) 
# This version of ellipse works, but you must set grad=F,is.there.hess=F 
# and plot.it=F 
# 
# Fit a least-squares ellipse centered at 0 with semi-axes (a, b) 
# and angle to the origin theta, to the data in x, y. The ellipse 
# is now parameterized by a, b (not e) and theta, in that order. 
# ais always > b; we can enforce that at the end. 
fr 

fet Semacrix(x) 6. ncol(x) > 1) { 

if (any (dimnames(x)[[2]] == "y")) { 
y <= x[, yt) 
if (any (dimnames(x)[[2]] == "x")) 


xX <- x; esc] 
else x <- x[, 1] 
} 
else { 
Veen ey eee 
ue ee | 
} 
} 
mm ( Ss. dist (x) ay 
if(any(names(x) == "y")) { 
y <- xS$y 
if (any (names(x) == "x")) 
xX <- xSx 
else x <- x[[1]}] 
} 
else { 
Ve flict) 
x <- x[[1]}] 
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} 


If a is supplied, and it's a vector, then we've been given 
starting points for all the parameters. Use ‘em, first making 
sure that there is the right number (3 if we're not fitting 
the center, and 5 if we are). These #’sS increase 1 for every delta 
estimated. 
if(!missing(a) && length(a) > 1) { 
if((fit.center && length(a) != 5) {| (!fit.center && 
length(a) != 3)) 
stop (paste("Parameter vector has length", length(a), 
", ©eXpecting “, i1felse(fit.center, 5, 3))) 
b <- a[2] : 
@o<— ss qr ia(b/a.(4)) *2) 
theta <- a[3] 
Li ( fat Center) if 
center.x <- a[4] 
center.y <- a[5] 


a <- afl] 
} 
else { 
if (missing(a)) 
ae<-w0eo eal ier (range <).) 
if (missing (b) ) 
Db (<=20n om “drt ii(cange (y)) 
lila aap) 
Switcheroo <- a 
a <- b 
b <- switcheroo 
} 
e <- sqrt(l - (b/a)%*2) 
if(missing(theta)) { 
Us. out@<- “I Stitt y) 
theta <- atan(ls.outScoef[2]) 
} 
center.x <- mean (xX) 
center.y <- mean(y) 
} 
tt <- ell.tt(x, y) 
ip lOr sr) a 
gtaph <—- dev.list() ["wimegraph™)] [1] 
if(length(graph) == || all(is.na(graph) ) ) 
win. graph () 
else dev.set (graph) 
} 
if (grad) 
Grad tung <— cil grad sir 
else gqrad.func <— NULLS 
if(fit.center) { 


start.vec <- c(a =a, b=b, theta = theta, center.x = 
CcCenter.x, center.y — center.y, delta = delta) 

Tower -Vece a— (G0 OOO 0d 0 2010 Oia em yee lke 
Lye 


WDPer. VEG = = ye (ene an) es ee eee ey oe ee) 
} 


else { 
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Suateen een ——etae— ia, > — Db, theta = theta, delta = delta) 


Owen. Vecw mc emuugdt, 0.00001, =2 * pl, - Inf) 
Uppetevec —— Fellini. Int, 2 * pl, Ins 
} 


out <— nilminb(start = start.vec, objective = ell.res.II, gradient 


= grad.func, hessian = is.there-.hess, lower 


| 


Lover Vee mmuUpPper = Wpper.vec, tt = tt, my.x = x, my.yY 
— SV, plotete — plot.1t, chat = chat, is.there.hess = 


is.there.hess, fit.center = fit.center, class.I 
class.I, which.type = which.type, step.min = 10 
-MachineSdouble.eps, scale.upd = 1) # 


+ ## p.-names <- names (out$parameters) 
### = cat("In ellipse, check p.names\n") 


— 
— 


O * 


# ## browser () 
1f(length (outSparameters) > 4) 
names (outSparameters) <- c("a", "b", "theta", "center.x", 
"“center.y", "delta") # Beats me... 
else names (outSparameters) <- c("a", "b", "theta", "delta") 


out$sigma.sq <- outS$obj/length (x) 
out$Ssigma <- sqrt(out$sigma.sq) 
outSaux <— NULL 
outSx <- x 
OUESY <- y 
if(length(outShessian) > 0) { 
1f(qr(outS$hessian)$rank < ncol(out$hessian) ) 
outScov <- "Can't invert Hessian" 
else outScov <- out$sigma.sq * solve(outShessian) 
} 
OUESTE <= EE 
SutSstitted- ce <— ell te(x = center.x, y -— center.y) 
pred<-ell.pred(out$Sfitted.tt, out$parameters["a"], 
outSparameters["b"],outS$parameters["theta"], 
return.unrotated.too = F, fit.center = fit.center, 


center.x = ifelse(fit.center, outSparameters["center. 
0), center.y = ifelse(fit.center, outSparameters["cent 
0)) 
outSfitted.x <- predSx 
outsSfitted.y <- preds$y 
Scot cued et e<— SGre(ipredox 2 + predSy~2) 
elass(out) <- "el#spse" 
Becurn (Oue) 
} 
> ellipse.III 
function(x, y, a = 2, b = 3, theta = 0, delta.a = 0, delta.b = 0, 
fit.center F, grad = T, is.there.hess = T, plot.it =T, c 


memeclass.£ = rep(lT, length(x)99 


This is for Class-III ellipses 
This version of ellipse works, but you must set grad=F,is.there 
and plot.it=F 


Fit a least-squares ellipse centered at 0 with semi-axes (a, b) 
and angle to the origin theta, to the data in x, y. The ellipse 
is now parameterized by a, b (not e) and theta, in that order. 
ais always > b; we can enforce that at the end. 
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sy] j 
ery, 


hat = 


-hess=F 


Lf (ts. Matrix (x )reee Neo mae > 1) a4 


if (any (dimnames (x) [[2]] == "y")) { 
4 <- x([, Ee 
if (any (dimnames (x) [(2]] == "x")) 
xX <- See ae? | 


else x <- x[, 1] 
} 
else { 
Van ly 2) 
x <- x[, 1] 
} 
} 


If (iscli2seiees 2 { 
if(any(names(x) == "y")) { 
y <- xSy 
if (any (names(x) == "x")) 


x <= 2S 
else x <- x[[1])] 
} 
else { 
Vee <1 2 1H) 
eee (h(i) | 


} 


e 
# If a is supplied, and it's a vector, then we've been given 
# starting points for all the parameters. Use 'em, first making 
# sure that there is the right number (3 if we're not fitting 
# the center, and 5 if we are). These #’s increase 1 for every delta 
# estimated 
if(!missing(a) && length(a) > 1) { 
1f((fit.center && length(a) != 5) || (!fit.center && 
length (a ) != 3)) 


stop (paste("“Parameter vector has length", length(a), 
", expecting ", ifelisie(fit.center, 5; 3))) 
b <- af[2] 
e <- sgqrt(i - (b/a[1])%*2) 
theta <- a[3] 
lf (£METecenter) { 
center.x <- a[4] 
center.y <- a[5] 
} 
a <- afi) 
} 
else { 
if (missing(a)) 
a <=) 0.5 ~Sdil fr range G.).) 
if (missing (b) ) 
be - 0-522 di fE (range (aay 
TE (a <b) { 
switcheroo <- a 
a <- b 
b <- switcheroo 
} 
e <- sqrt(l - (b/a)%*2) 
if(missing(theta)) { 
TS. OUE =~ Sls fie ow ey) 
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# # # 
HHT 
# # # 


# HH 
HHH 
# # # 


theta <- atan(ls.outScoef([2] ) 
} 


Gemeer 6 <— Meant x) 
€etieer.y <- mean (y) 


fee Cll Ge x, yy) 


PaplouaLte) -{ 
Geapn <— dev ist () ["win.graph"™] [1] 
if(length(graph) == Pitail (is na(gqraph))) 


win. graph () 
else dev.set (graph) 
} 
if (grad) 
grad.func <- ell.grad.III 
else grad.func <- NULL 
Pi(tit. center) 4 
start.vec <- c(a =a, b= b, theta = theta, center.x = 
center.x, center.y = center.y, delta.a = delta.a, 
delta.b = delta.b) 
hOWe be ema en O  CUOOT SO vC000 lye 2 a-spi —sinnt, | -— Int, 0; 
0) 
UpDeGwiecm<— CULNE, SINt,. co “Dismint, Inf, int, Int) 
} 
else { 
Start.vec <- c(a =a, b= b, theta = theta, delta.a = 
delta.a,delta.b = delta.b) 
teweoevee.<-— c(0.00001, 0200001, -2 * pl, 070) 
Uppewe ween a C(Int, inf, 2 * pi, Inf, Int) 
} 
out <- nlminb(start = start.vec, objective = ell.res.III, gradient 
= grad.func, hessian = is.there.hess, lower = lower.vec, 


upper = upper.vec, tt = tt, my.x = x, my.y = y, plot.it = 
plot.it, chat = chat, is.there.hess = is.there.hess, 
fit.center = fit.center, class.I = class.I, step.min = 100 * 


.Machine$double.eps, scale.-upd = 1) # 
p.names <- names (outSparameters) 
cat("In ellipse, check p.names\n") 
browser () 
if(length (outSparameters) > 4) 
names (outSparameters) <- c("a", "b", "theta", "center.x", 


"eenter.y", "delta.a", “delttar.b”) # Beats me... 
else names(outSparameters) <- c("a", "b", "theta", "delta.a", 
"delta.b”) 


outS$sigma.sq <- outS$obj/length (x) 
outS$sigma <- sqrt(outSsigma.sq) 
outSaux <- NULL 
SutSx <= x 
SGUcoy <—)¥ 
if(length(outShessian) > 0) { 
if(qr(out$hessian)$rank < ncol(out$Shessian) ) 
outScov <- "Can't invert Hessian" 
else outScov <- outSsigma.sgq * solve(outShessian) 
} 
SUtote <- Ef 
outSfitted.tt <- ell.tt(x - center.x, y - center.y) a 
pred <- ell.pred(outS$fitted.tt, outS$parameters["a"], 
outSparameters [ 
"b'], outSparameters["theta"], return.unrotated.too = F, 


ay 


# # # fit.center = fittzvcenter, center.x = ifelse(fit.center, cw 
# # # parameters["center.x"], 0), center.y = ifelse(fit.center, 
+ # # outSparameters["center.y"], 0)) 
# + + outSfitted.x <- pred$x 
# # + out$fitted.y <- predSy 
# # # Outstitted 9 £.<— sare (oredox 2 + predsy-z) 
class(out) <- "ellipse" 
return (out) 


> ell.res 
function (params, Ct,7 my.x, my.y, 1s.there Ness, fit.center, plct 2 
chat) 


{ 
# 
# ell.res: Compute objective to be minimized. 
4 
# This computes the objective function: the sum of squared 
# differences between the observed points on the ellipse 
# (after transformation) and the predicted ones. 
# 
# "params" is the vector (a, e, theta). Get them out, and 
# compute rat, the ratio a/b. 
¥ 
a <- params[1] 
e <- params[2] 
shag (Ses ae SIe)) 
return (1000) 
be <= a, ~*~ SGemite-— e222) 
theta <- params[3] # 
if(fit.center == T) { 
center.x <- params[4] 
center.y <-— params[5] # 
CU <—-<eChi ery 2s = Senter. xX,2my-) = center.) 
} 
else { 
CGnbers wea, COhnecr.y <- 0 
} 
fitted.r <- ell pred(tt, a, b, theta, £10.center = £1 center, 
center.x = center.x, center.y = center.y) # 
new.x <- fitted.r$x 
new.y <= fitted. xrs5y 4 
# Plot it: add dotted lines at x = 0 and y = O. 
i 


Peto LoOt. 1 ti) 4 
plot (my.x, my.y, Xlim = range(my.x, new.x), ylim = 
Fange iy ele 
abline(h = 0, l1ty = 2) 


abline(v = 0, lty = 2) # 
points(new.x, new.y, pch = 1, col = 4) # 
poin@sH center .x, center. y, pel — |, col = 2) 


} 
+ # # for(i in i:lengeiee.) 


# # # polar(ele7 O45), repice la 2) se evpes— es uls) # 
A 

# Get fitted x and y; compute and return objective. 

5 
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Obi eounny <a etewe 92) + Ssum((My.y - new.y)*2)  # 


PCC at oa teed wc mgmeriicd, 4),) ',b:", signif(b, 4), “",th:", 
SrOtert(Ehetamme mere Se(f1t.center, paste("7x,y:", signif ( 
Cenvere, 6°) Gest ammrercenter.y, 4)), ""), “Fob j:", 
SuGmtf(Ob) 2 2) wenn) ## # cat ("BIWee") 


# # # ell.grad(params, tt, my.x, my.y, is.there.hess) 
mee |6CCAE("\n") 
return (obj) 


> ell.res.II 


Mmrietion(params, tt, My.x, My.y, is.there.hess, fit.center, plot.it, 


GhatewclaSs s,s welch .type) 
ell.res: Compute objective to be minimized. This version is the 
wirass-Ii one. 
This computes the objective function: the sum of squared 
Mireerenices DeEween tne @bserved points on the ellipse 


(after transformation) and the predicted ones. 


"params" is the vector (a, b, theta). 


$6 SE OSE SF SE OSE OSE OE HEHE 


a <- params[1] 

b <- params[2] 

theta <- params [3] 4 

mnt is. Wwe Lass. 1) ) -{ 
class.I <- rep(T, length (my.x) ) 
delta <- 0 

} 


if(fit.center == T) { 
center.x <- params [4] 
center.y <- params[5] # 


delta <- params [6] 

Cie eclMmNcCNMvVene- CencCer.x,; My.-y -— center.y) 
} 
else { 

center.x <- center.y <- 0 

delta <- params [4] 


} 
mExSum(class. 1jye< length(my.x)) { 


PLGread- 6-1. <- ell.pred(tt(class.1), a, Db, theta, -f£1t-center 


= fit.center, center.x = center.x, center.y = 
center.y) # 

if(which.type == 1) 
EVecederaiie<— ell pred(tt,! class .i)7.a.7 aelea, 


theta, fit.center = fit.center, center.x = 
Center.xX, Center.y = center.y} 

else fitted.r.II <- ell.pred(tt[!class.I], a, b + delta, 
theta, fit.center = fit.center, center.x = 
center.x, center.y = center.y) # 


} 
else) £fitted-r.1 <= ell preditt, a, b, theta, fic.center «= 


fit.center, center.x = center.x, center.y = center.y) # 
new.x <- numeric(length(my.y) ) 
new.y <- numeric(length(my.y) ) 
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new.* (class 20 lief) ttecusa io 

new.y (class. Wie. at teed. doy 

ifisum(!class 1 2 10). ( 
new.x{!{class.ii))<— fitted.r.L1isx 
new.y[!class.I]) <- fitted.r.IIsSy 


ft 
fr 
# If plot it, add dotted lines at x = 0 and y = O, plus points. 
fr 
iGO Gee) aed 
pPLOt(my-x,emy.y, xlim = range (my.x, new.x), ylim — 
range (my-y, new. y) ) 


abline(h = QO, l1ty = 2) 
abline(v = 0, 1lty = 2) # 
points (new.x, new.y, pch = 1, col = 4) # 
points (center.x, center.y, pch = 1, col = 2) 
} 
it 
Gat {("grad.norm as 9", sum(ell.grad.liparams, “ct, *my.x, my*y, 
is.there.hess, fit.center, class.I, which.type)%*2), "\n"™) 
4 
# Get fitted x and y; compute and return objective. 
# 
Ob] <- Sum((my.x - new.x)*2) + sum((my.y —- new.y)~Z) # 
LeEtochnats 


cat("a:", signif(a, 4), ", delta: ", signif(delta, 4), 
“7b: ",Signititb, 4), “,th:", signif(tneta, 4). 1felsee 
fit.center, paste(";x,y:", signif(center.x, 4), 
Sagmin(center:y, 4)), ""“), "70b3:", Signak(obj], an 
WN) 


return (obj) 


> ell.res.III 
function (params, Et, my.x, My-y, is.there-hess, fLit.center, plotwase 
chat, class.TI) 


ell.res: Compute objective to be minimized. This version is the 
Class-III one. 

This computes the objective function: the sum of squared 
differences between the observed points on the ellipse 


(after transformation) and the predicted ones. 


"params" is the vector (a, b, theta). 


+e SE Se OEE HEHEHE 


a <- params[1] 

b <- params [2] 

theta <- params [3] # 

LE(1iSsnul (clase ae 
class.i <— rep (Tl ,tengtn (myx) ) 
delta.a <- delta.b <- 0 


} 
if(fit.center == T) { 
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center.x <- params[4] 

center.y <- params[5] # 

delta.a <= params[6] 

delta.b <- params[7] 

Cee crete Wye n> Center.x, my.y - centers 
} 
else { 

Genter,<.<- Center.y <—- 0 

delta.a <- params [4] 

delta.b <- params[5] 


} 


if(sum(class.I) < length(my.x)) { 
Pimecamnret.=-— cll. pred(ttiiclass.I], a, b, theta, fit.center 
= fit.center, center.x = center.x, center.y = 
center.y) # 


Piece liv ell,prea(cr(!class.I)], a + deélta.a, b + 
delta.b, theta, fit.center = fit.center, center.x 
center.x, center.y = center.y) # 

} 
else fitted.r.I <- ell.pred(tt, a, b, theta, fit.center = 
fit-center, center.x = center.x~ecenter.y = center.y) 
new.xX <= numeric(length (my. y) ) 
new.y <- numeric(length(my.y) ) 
new.x[class.I) <- fitted.r.I$x 
new.y(class.I] <- fitted.r.ISy 
Pecsum(!class~1) 2-0) 
new.x[!class.I] <- fitted.r.IIS$x 
new-yi'teclass.1]' <— fitted.r. Lisy 


4 
f 
# If plot it, add dotted lines at x = 0 and y = 0, plus points. 
# 
et (pLOeeae) 4 
plot (my.x, my.y, xlim = range(my.x, new.x), ylim = 
Fancge(Nyny,, mew. y)) 
abline(h = 0, l1ty = 2) 
abline(v = 0, lty = 2) # 
points (new.x, new.y, pch = 1, col = 4) # 
Poumes(centen. weeenter.y,-pch = 17 col = 2) 
} 
+ # # 
## Gat Grad snorm is ", sum(ell.grad.11 (params, tt, my.x, my:y, 
# # is.there.hess, fit.center, class.I, which.type)*2), "\n") 
+ 


# Get fitted x and y; compute and return objective. 


obj <- sum((my.x - new.x)*2) + sum((my.y - new.y)*2) # 
mea Chat) 

Gat( “an, sSsignitvta, 4), ", delta.a: Waesignit (delta.a, 
Seem, SIiGhLE(b, 4)eeeidelta.b: ay signif (delta.b, 
 ,en: -» Signit(theta, 4); ifetse(fit.center, paste ( 
oy 7 SLOWMLE (Center.x, 4), osigniésé (center. ypu) ), 
wee 2 OD): 7 -SLONL SO), 4), \n- 


return(obj) 
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> ell.grad 
function(params, tt, my.x, my.y, is.there.hess, f£1t.center) 


{ 


HF 
# HH 
HHH 


tol <— le-006 
a <- params([1] 
e <- params([2] 
theta <- params [3] 
if(fit.center == T) { 
center.x <- params(4] 
center.y <- params([5] 
bbc ell be my. Xo Cente tx seiny ec CNeECIn wy) 
} 
else { 
Center.x <—- céenter.y <= 0 
} 
b <—"a = Sqrb Gre eZ) 
fitted <— ell predict, a, b, theta, Lreturn.unrotatea.too = T, 
fit.center = fit.center, center.x = center.x, center.y = 
Ccenteue) 
xprime <-— fitted$x.prime 
yprime <- fittedSy.prime 
x <—- fittedsx 
y <- fittedSy 
cos.theta <- cos(theta) 
Sin.theta <- sin(theta) 
Cos. Z2.tt.theta <— cos(Z * (tt. — tEheta)) 
Sime2.tt. theta <— sin(2Z)~ Sct — theta) } 
Sinsg.tt.theta <- (sin(tt -— theta) )%*2 
cossq.tt.theta <- (cos(tt - theta))%*2 
Sinsq.2.tt.theta <— “(sin(Z2 * (tt - theta)))“*2 
Consd.2.ttstmetal <=—(cos(2 ~~ (Et — theta)jy4Z 
one Minus). sqy—— i —- e€%2 
denom <- cossq.tt.theta * one.minus.e.sq + sinsgq.tt.theta 
dxprime.da <- xprime/a 
dxprime.de <- (( - a*2/(4 * xprime)) * (e * 
sinsq.2.tt.theta) )/ (denom* 2) 
dxprime.de[abs(xprime) < tol] <— 0 
dyprime.da <- yprime/a 
dyprime.de <- - (one.minus.e.sq * xprime * dxprime.de + e * (a%2 
-xprime%*2) )/yprime 
dyprime.de[abs(yprime) < tol] <- 0 
dx day<— cos theta“ tadxprimesda — sin-eneta ~*~ adyprime sda 
dx.de <- cos.theta * dxprime.de - sin.theta * dyprime.de 
dyadaG-= Simmachetas, sOxprame da +tryecos. theta SeaVpraime-da 
dy.de <- sin.theta * dxprime.de + cos.theta * dyprime.de 
x.diffi <- my.x - x 
y.diff <- my.y - y # 
grad.mat <- matrix (0, Vengths).. 3) 
grad .matly) 1) == 2. (x dare Cx dat yaar tl) “ayaa 
Grad. Ma tie le oe lI ee Cle Glew cll bien acy dc) 
GrYadsa <Seo2ee sum (xX. atit dx. daw -cdltif. —dy7da) 
Grad. @ =a Sum (x. GLEE * dx-de = ect ft = adv de) 


Himes One wMatnusme Sos sisin(2 ~ te —eeencta) ) 

dxprime.dtheta <= (a"°2/(2 * xprime)) *  (num/denem 2) 
dxprime.dthetafabs(xprime) < tol] <- 0 

dyprime.dtheta <- - (one.minus.e.sq * xprime * dxprime.dtheta) / 


yprime 
dyprime.dtheta[abs(yprime) < tol] <- 0 
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Gxea cet. ei mee cemiatm 7) + cos.theta * dxprime.dtheta - 
Sin.theta * dyprime.dtheta 
Oy -Gthetame §(X)—scenter,x) +isin.theta * dxpramemdtheta + 
cos.theta * dyprime.dtheta 
grad-theta <=—-2°* Sum(xmeiftft * dx.dtheta + y.diff * dy.dtheta) 
# 
# # # Glad tcc Aeeemairt ~sdxedtheta +- y.difft * dy.dtheta) 
# # # Gat ("Grad mat approx. 152. -\n") 
## # print(t(grad.mat) %*% grad.mat) 
if(fit.center == F) 
grad <— c(guacma, grad-e, grad.theta) 


else { 
dxprime.dt <- - dxprime.dtheta 
dyprime.dt <- - dyprime.dtheta 
RoC my sie mCeheCm a a7. + (myoy — center ly) 2 
dit a0 S<= my yemwconter.V)/R.s¢ 
dt.dyO <- - (my.x - center.x)/R.sq 
dxprime.dx0O <- dxprime.dt * dt.dx0 
dxprime.dyO <- dxprime.dt * dt.dy0 
dyprime.dx0 <- dyprime.dt * dt.dx0 
dyprime.dyO <- dyprime.dt * dt.dyO # 
dyprime.dxprime <- one.minus.e.sq * ( - xprime/yprime) 
dyprime.dxprime[abs(yprime) < tol] <- 0 
ax.dx0 <— (cos.theta * dxprime.dx0) — (sin.theta * 
dyprime.dx0) + 1 
dy dx0s<-s1sino-theta * dxuprime.dx0)) + (cos-theta * 
dyprime.dx0) 
dx.dy0Q <- (cos.theta * dxprime.dy0O) - (sin.theta * 
dyprime.dy0) 
dy.dyO <- (sin.theta * dxprime.dy0) + (cos.theta * 
dyprime.dyO) + 1 
Guacd. x<Omia—=—2e" SUMN(x GLEtT = dxsdx0 7 avec fl “dy dx) 
GeaceyUh- eer vSuntx, Grin GxsavOet y.OLifetsdy.dyo} 
Grad <— c(grad-a, gradier grad.theta, grad.x0, grad.y0) 
} 
if(is.there.hess == F) 


return (grad) 
d2xprime.da2 <- d2yprime.da2 <- 0 
d2xprime.dade <- dxprime.de/a 
d2yprime.dade <- dyprime.de/a 
d2xprime.dadtheta <- dxprime.dtheta/a 
d2yprime.dadtheta <- dyprime.dtheta/a 
ddenom.de <- -2 * e * cossg.tt.theta 
ddenom.dtheta <- - e*%2 * sin(2 * (tt - theta) ) 
terml,<- ( - a*2 * sinsg.2.tt.theta) /4 
xprime.denom.sq <- xprime * denom*2 
ix prRiIMendey <—- xpreime.aGeiom.sgq —- cc ~ (axprime.de “~saenom 2 + 2° 
xprime * denom * ddenom.de) 


d2xprime.de2 <- (terml/xprime.denom.sq*2) * d2xprime.de2 

dZ2xprime.deZ{abs(xprime) < le-0Q06] <- 0 

d2yprime.de2 <- one.minus.e.sq * (-2 * xprime * d2xprime.de2 - 2 * 
dxprimerde 2). +76 *.e * xprame * daomrme.de — 25%=(a,2°— 
xprime%*2 ) 

Germl <= .— a 25" e 

d2xprime.dedtheta <- - xprime.denom.sq * sin(4 * (tt - theta)) - 


Sinsqd,2.ce,eheta * (xprime ~*~ denom * daadenom.dtheta + 
dxprime.dtheta * denom%’2) 
d2xprime.dedtheta <- (( - a*2 * e)/xprime.denom.sq%*2) * 
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d2xprime.dedtheta 
Geese ime. Gedcenera[abs(xprime j= 1e-U0G).— 0 
d2yprime.dedtheta <- (-l/yprime*2) * (yprime * ((one.minus.e.sq 4 


xprime * d2xprime.dedtheta + dxprime.dtheta * dxprime.de)) - 
2 * e * xprime * dxprime.dtheta) - one.minus.e.sq * xprime * 
dxprime.dtheta * dxprime.de) 

d2yprime.dedtheta[abs(yprime) < le-006] <- 0 

dnum.dtheta <- —2Z “sone minus.e€:-sq-— "cos (2 * (te —-seneta 

d2xprime.dtheta2 <- xprime.denom.sq * dnum.dtheta - num * (2 * 
xprime * denom * ddenom.dtheta + denom*2 * dxprime.dtheta) 

d2xprime.dtheta2 <- (d2xprime.dtheta2 * a%*2)/(2 * 
xprime.denom.sq*2)d2xprime.dtheta2[abs(xprime) < le-006] <-0O 

d2yprime.dtheta2 <- - (one.minus.e.sq/yprime) * (yprime * (xprime 
* d2xprime.dtheta2 + dxprime.dtheta%*Z) - dyprime.dtheta * ( 
xprime * dxprime.dtheta) ) 

d2yprime.dtheta2[abs(yprime) < le-006] <- 0 

d2x.dade <- cos.theta * d2xprime.dade - sin.theta * d2yprime.dade 

d2y.dade <- sin.theta * d2xprime.dade + cos.theta * d2yprime.dade 

d2x.de2 <- cos.theta * d2xprime.de2 - sin.theta * d2yprime.de2 

d2Zy.de2 <= Sin.theea “ d2xprime.deZ + cos.theta * dZyprime.dez 

gGrad.aZ2 <= 2 * sum(dx.da%Z + dy.da%2) 

Grad.ae <- 2 * sum( —- x.diffi * d2x.dade + dx.da * dx-.de - y.dijG@ime 
dZ2y.dade + dy.da * dy.de) 

d2x.dadtheta <- cos.theta * d2xprime.dadtheta - sin.theta * 
d2yprime.dadtheta - dy.da 

d2y.dadtheta <- sin.theta * d2xprime.dadtheta + cos.theta * 
d2yprime.dadtheta + dx.da 

d2x.dedtheta <- cos.theta * d2xprime.dedtheta ~- sin.theta * 
d2yprime.dedtheta - dy.de 

d2y.dedtheta <- sin.theta * d2xprime.dedtheta + cos.theta * 
d2yprime.dedtheta + dx.de 

d2x.dtheta2 <- cos.theta * d2xprime.dtheta2 - sin.theta * 
dZyprime.dtitetaZz —~2 * dy.dtheta @ x 

d2y.dtheta2 <- sin.theta * d2xprime.dtheta2 + cos.theta * 
dZyprame=dtheta2 + 2 *7dx.dtheta +7 y 

grad.atheta <- 2 * sum( - x.diff * d2x.dadtheta + dx.da * 

dx.dtheta —- y.diféi * d2y.dadthetarmee dyeda * dy.dtheta) 

GuadeeZ (<— 927 * sume «.diff * (dZxsde2u eax dec2. yd btu 
d2y-de2 + dy-de2} 


grad.etheta <- 2 * sum( - x.diff * d2x.dedtheta + dx.de * 
dx.dtheta - y.diff * d2y.dedtheta + dy.de * dy.dtheta) 
Grad: thetaZ <= 2° * =sum(e—sx Gilt tuedZx dthetazet dx. dtncta 2a. 


y-diff * d2y.dtheta2 + dy.dtheta’*2) # 
tt grad.mat <- matrix(c(grad.a2, grad.ae, grad.atheta, grad.ae, 
gGrad.e2, grad.ethetay, qradeatheta, gnad.ctheta, 
grad.thetaZ), 3,7sen) 


Hitt Gabi ressil ane dep rox. sls ee) 
+ # # print (solve(grad.mat) ) 
bf (ie secnter => Fb) 


hessian <- c(grad.a2, grad.ae, grad.e2, grad.atheta, 
grad.etheta, grad.theta2) 
else { 


Second derivatives: a and x0, a and y0 


tHe sth te Sk 


d2xprime.dadx0 <- dxprime.dx0/a 
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d2yprime.dadx0 <- dyprime.dx0O/a 

d2xprime.dady0O <- dxprime.dy0/a 

d2yprime.dady0O <- dyprime.dy0/a 

d2Zx.dadx0@-emcos, theta * d2xprime.dadx0 - sin.theta 
d2yprime.dadx0 

d2y.dadx0 <- sin.theta * d2xprime.dadx0 + cos.theta 
d2yprime .dadx0 # 


# and the gradient 


+e te 


Grad axU0ma a2 wsum(=-— xX2-diff * d2x.dadx0 + dxida * 
Vaewer = G2y dadx0 - dy dame ay-ax0) 

d2x.dadyO <- cos.theta * d2xprime.dady0O - sin.theta 
d2yprime.dady0 

d2y.dady0 <- sin.theta * d2xprime.dady0 + cos.theta 
d2yprime.dady0 7 


Giada omnes 2a Sun (= x.dLfi *ed2x.dady0 + dx.da * 
Wachee G27 dadyU-+ dy.da * dy.dy0) + 


ddenomaegx0 <- e“Z * sin.2.tt.theta * dt.dx0 
ddenoimacy0 <— 6€°2 * Sin.2.tt.theta * drsdy0 
A <- xprime * denom*2 


@A.dx0 <- denom * (2 * xprime * ddenom.dx0 + denom * 


dxprime.dx0) 


qadA.dyO <- denom * (2 * xprime * ddenom.dy0 + denom * 


dxprime.dy0) 
out.front <- - (a*2 * e)/4 
z <- sinsg.2.tt.theta/A # 
z[abs(A) < tol] <- 0 


* 


* 


ax.dax — 


* 


Gil. 0 a 


2p eine -Cedmaeoutl.f rome ~ ((2°* sin(4 * “tt — theta)) * 


at .adx0)/A— ((dawax0 * z) /A)) 


GO2xprime-~cdeady0. <— GUutL.Tromge.- ((2)* Sin(4 estes theta) ) * 


dt.dy0)/A - ((dA.dyO * z)/A)) 
d2xprime.dedxO[abs(A) < tol] <- 0 
d2xprime.dedyO[abs(A) < tol] <- 0 # 


# Here's one from Mathematica. 


# 


# 


d2yprime.dedx0 <- ( - (3 * e)/2 * yprime * dt.dx0O * 
sin.2.tt.theta) /denom*2 

d2yprime.dedx0O[abs(yprime) < tol] <- 0 

d2x.dedx0 <- cos.theta * d2xprime.dedx0 - sin.theta 
d2yprime.dedx0 

d2y.dedx0 <- sin.theta * d2xprime.dedx0 + cos.theta 
e2yprime.dedx0 

gGrad.ex0 <— 2 “ssum( — x ditt * d2x,dedx0 4 .dx.de * 
YVeamnt ~*~ dZy.dedx0 + dy.de * dy.dx0) 

d2yprime.dedyO <- ( - (3 * e)/2 * yprime * dt.dy0O * 
sin.2.tt.theta) /denom*2 

d2yprime.dedyO[abs(yprime) < tol] <- 0 

d2x.dedy0 <- cos.theta * d2xprime.dedy0O - sin.theta 
d2yprime.dedy0 

d2y.dedy0O <- sin.theta * d2xprime.dedy0O + cos.theta 
d2yprime.dedy0 


* 


* 


dx.dx0O - 


(tad ,6 Uso SU owee ttf ~a2Ze, dedy0. tod de, 80x. dy07— 


y.diff * d2y.dedy0O + dy.de * dy.dy0) 7 


# Here's another from Mathematica. 


# 
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OU. ErOnE <= Sizprines= sii —92 * e° 254 €°2) 
€0s5.2-Et~ thera), denomeZ 

out.front[abs(denom) < tol] <- 0 

d2xprime.dthetadx0 <- out.front * dt.dx0 

d2xprime.dthetady0O <- out.front * dt.dyO # 

out.front <- ((one.minus.e.sq) * yprime * (3 - 2 * denom))/ 
denom”*2 

out.frontf{abs(denom) < tol] <- 0 

adZ2yprime-denetadxO.-- out.front * dt.dxo 

d2yprime.dthetadyO <- out.front * dt.dyO # 


d2x.dthetadx0 <- - dy.dx0 + cos.theta * d2xprime.dthetadx0 
- sin.theta * d2yprime.dthetadx0 
GQ2ye dtnecadz0e<— (dx dx0 —- 1)e+ sin. thera = 
d2xprime.dthetadx0 + cos.theta * d2yprime.dthetadx0 
gnad-thetax0)<— 2 5* )siim( — x2difti * d2x-dthetadx0 + 
deeadtheta * "dx: -axOg—ay.dittiv-sdZy.dthetadx0 4 
dy -dtheta *sdy2d~0) 
d2x.dthetadyO <- - (dy.dy0O - 1) + cos.theta * 
d2xprime.dthetady0O - sin.theta * d2yprime.dthetady0 
d2y.dthetadyO <- dx.dy0O + sin.theta * d2xprime.dthetady0O + 
cos.theta * d2yprime.dthetady0 
Grads thetayO.<—- 2 * ‘sum( — xX.dil£i ~ dz2x-dthetadyo0 + 
dwraemetra dxcay0 — yodiff * dZy.dthetady® + 
ay -dtneta * dy-dyo0) 
E 


Zia as OZ eae Zeenat ~cir0 7) diy dy!) 

dZbax0dy Oi — ey RSG etree Gt ax 0) Z 

d2t.dy02 <- - d2t.dx02 

aZxprime. axOZs<— ) - cxprimne-arneta - dzt-ax0Z — 
d2xprime.dthetadx0 * dt.dx0 

adZyprime.ax02 <—) — dyprime-dencea ~saZt-dx0Z — 
d2yvyprime.dehnetadx0) * deed 

d2x-ax02 -<— cos.theta,~ dzZxprime.dx@z — sin. theta * 
d2yprime.dx02 

ad2y.ax0Z2 <- = sin.theta * dZxprime-dx02-+ cos theta * 
d2yprime.dx02 

grad.x02Z2 <= 2 * sum( =e -diff * d2x7dx02 4+ °dx.dx0°2 = y 7disme 
* GZ. ax02 = tay = cin): Zum 


d2xprime.dxOdy0O <- - dxprime.dtheta * d2t.dx0Ody0O - 
d2xprime.dthetady0O * dt.dx0 
d2yprime.dxOdyO <- - dyprime.dtheta * d2t.dx0dy0 - 


d2yprime.dthetady0O * dt.dx0 

d2x.dxO0dy0 <- cos.theta * d2xprime.dx0Ody0 - sin.theta * 
d2yprime.dx0dy0 

dZy.dxOdyO <= sin.theta * dZxprime.dx0dy0O + cos.theta * 
d2yprime.dx0dy0 

Grad.x0vyO <—92 2 sum (@—= = dit ie eczx axle Ont dx.ax0 * 
dx. dy@™—- y.diftf *"d2y.dx0aye J say-an0e. dy ay0) ¢ 


d2xprime.dy02 <- - dxprime.dtheta * d2t.dy02 - 
eZxprime athnetady?-* denayo0 
dZyprime.dy02 <= — dyprime.dtheta = tered, 02 — 


d2yprime.dthetadyO * dt.dy0 
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tH # 


G2x.aVvUce—-wees-tneta * d2Zxprime.dy0OZ -— sin.theta * 
d2yprime.dy02 
a2yedy02e<— Siac heta * dZxprime.dyO2m+ cos.theta * 
d2yprime .dy02 
Geaa=yO2 <-scmmeecuni es xsditt * dzZx.dyO2 + Gx.dy0*2 - y.diff 
Sea2V-cdvyuz-+ dy.ayo°2Z} 
hessian <-— c(grad.a2, grad.ae, grad.e2, grad.atheta, 
grad.etheta, grad.theta2Z, grad.ax0, grad.ex0O, 
grad.thetax0, grad.x02, grad.ay0O, grad.ey0O, 
grad.thetay0O, grad.x0y0, grad.y02) # 
print (hessian) 
} 
thing <- list(gradient = grad, hessian = hessian) 
return (thing) 


> ell.grad.ItI 
GamectciLon (params, tt, my.x, my-y, is.there.hess, fit.center, class.I, 


{ 


which.type) 


tol <- le-006 a 
a <- params[1] 
b <- params [2] 
e <- sqrt(l - (b/a)%*2) 
theta <-— params[3] 
if(fit.center == T) { 
center.x <- params [4] 
center.y <- params[5] 
delta <-— params[6] 
Cees e Ei Cem ex Center x, my-ye- center.y) 
} 
else { 
delta <- params[4] 
GEenter.xe-— Centcr.y <- 0 


} 


if(sum(class.I) == length(my.x)) { 
fitted <- ell.pred(tt, a, b, theta, return.unrotated.too = 
T, fit.center = fit.center, center.x = center.x, center.y = 
Center, 


xpxrime <— fitted$x-.prime 
yprime <-— fitted$Sy.prime 
x <- fitted$x 
y <- fittedSy 

} 


else { 
flected L <— 7clleored(ttl[elass.[)], a, b,wtheta, 
return.unrotated.too = T, fit.center = fit.center, 
center.x = center.x, center.y = center.y) 
if(which.type == 1) 


fieced.II <—- ell. prea(tte|!class.I], a + delta7ab, 
theta, return.unrotated.too = T, fit.center 
= fit.center, center.x = center.x, center.y 
= center.y) 

else fitted.II <- ell.pred(tt[!class.I], a, b + delta, 

theta, 
return.unrotated.too = T, fit.center = 
fit.center, center.x = center.x, center.y = 


6] 


# 


# Okay. Here's 


# 


center.y) 


xprime <- yprime <- x <- y <- numeric(length(my.x) ) 


xprime[class.I] 


<= E£itted-ISx. prime 


xprime[!class.I] <- fitted.II$x.prime 
yprime{class.I] <- fitted.I$Sy.prime 
Vorime | velass. time fitted 1isy. prime 
x{class.I] <- fitted.I$x 
x{!class.I] <- fitted.II$x 
y{class.I] <- fitted.ISy 
y{!class.I] <- fitted.IISy 


} 


if(which.type == 1) { 


ana Geo (ay 


length (my.x) ) 


af{!class.I] <- a[!class.I] + delta 


} 


else { 


b <~ rep(b, length(my.x) ) 
b[!class.I] <- b[{!class.I] + delta 


} 


cos.theta <- cos(theta) 
Sin.theta <- sin(theta) 


cos.2 
Sine. 


Sinsq.tt.theta <- 
COSSQG- Ge eieea— 


-tt.theta <- cos(2 * 
tt.theta <-— sin(2 * 


Sinsd.2 tt. theta <— ~(Sin(Z * 
consg.Z.tt.theta <= (cos(2 * 
one.minus.e.sq <- 1 - e%2 


(tt - theta) ) 

(tt -— theta) ) 
(sin(tt - theta))%*2 

(cos(tt - theta))%*2 

(Et ={theta) yz 
(tt - theta)))%2 


denom <- cossq.tt.theta * one.minus.e.sg + sinsgq.tt.theta 


dxprime. 
dxprime. 
sinsq.2Z.- 
dxprime. 
dyprime. 
dyprime. 


dyprime. 


da <- xprime/a 
de <- 
tt.theta) )/ (denom*2) 
defabs(xprime) < tol] <- 
da <- yprime/a 

de <- - (one.minus.e.sq 
xprime’*2) ) /yprime 
defabs(yprime) < tol] <- 


where we go from "e" to 


(( - a*2/(4 * xprime)) * 


de. 
cox 
Gx 
aise: 


dy. 
ay: 
ay. 


db 
da 
de 
db 


da 
de 
db 


<— 


de. 


g<- 
g~— 
<- 


de. 
X.adiff <- 


- b 
cos 
cos. 
cos. 
db 
Sin. 
Sin. 
Sine 
db 

my - 


y-diff <- my. 
Gradza- andy 


grad.a <- sum(grad. 


# HH 


gradib <- sumi(Gnuace 


fiace 


.theta 


theta 
theta 


theta 
theta 
theta 


oe 
Vee ey, 
B= 8 


grad.e <- -2 
Glad. band 7 oe 


if(which.type == 
grad.delta <= sum(grad.a.indiv{[!'class.1I]) 


* 
* 
* 
* 


a 
* 
x 
b 


) 


Si 


dxprime. 
-de 
.de 


dxprime 
dxprime 


dxprime. 
. de 
.de 


dxprime 
dxprime 


(oxeclaate te ack 
; Tea) 
Sumtx, a1 ff a 
(ceca fi Faicisx 
nie) 


da 


da 


62 


ea 


0 


* xprEame ~ Gxprime.de + e * 


O # 


"wh ve 2 


Sin.theta * 
Sin.theta * 
de.db - sin. 


cos.theta * 
cos.theta * 
de.db + cos. 


-da + yeewitt * dy.da) 


dyprime. 
dyprime. 


theta * 


dyprime. 
dyprime. 


theta * 


(a%2 


We still need dx/y.de. 


da # 

de 
dyprime.dems 
da # 

de 


dyprime.de * 


Gx-Gew- VY -Giree Gy. ce) 


sabi + y ocdat f yen, di> ) 


else grad.delta <- sum(grad.b.indiv[!class.I]) 

num <- one.minus.e.sq * sin(2 * (tt - theta)) 
dxprime.dtheta <- (a*%2/(2 * xprime)) * (num/denom“%2) 
dxprime.dtheta[abs(xprime) < tol] <- 0 


dyprime.dtheta <- (one.minus.e.sq 
yprime 
dyprime.dtheta[abs(yprime) < tol] <- 
dx.dtheta <- (y — Ccenter.y) + cos 
Sin.theta * dyprime.dtheta 
dy.dtheta <- (x - center.x) 


cos.theta * dyprime.dtheta 


* xprime * dx¥prime atheta) / 


0 


-theta * dxprime.dtheta - 


+ sin.theta * dxprime.dtheta + 


grad thnebame.)—2 .*esum(xX-ditf * dx.dtheta + y.dififie* dy.dtheta) 


frame 


# Pein exists (Kil ilme", 
browser () 


if (fit.center 


== ) 


0) ) 


grad <- c(grad.a, grad.b, grad.theta, grad.delta) 
else { 

dxprime.dt <- - dxprime.dtheta 

dypamme.dt <- — dyprime.dtheta 

Reed amy. — center.) "2 + (my.y = center.y)°2 

dt.dx0 <- (my.y - center.y)/R.sq 

de-dy0"<—- ~—' (my.x = center 2a%/R.sq 

dxprime.dx0 <- dxprime.dt * dt.dx0 

depeime .dyO <—- dxprime.dt * dt .dyv0 

dyprime.dx0 <- dyprime.dt * dt.dx0 

dyprime.dy0 <- dyprime.dt * dt.dyO # 

dyprime.dxprime <- one.minus.e.sq * ( - xprime/yprime) 

dyprime.dxprime[abs(yprime) < tol] <- 0 

dx.dx0 <- (cos.theta * dxprime.dx0) - (sin.theta * 
dvorime cas) ae 

dy.dx0 <- (sin.theta * dxprime.dx0) + (cos.theta * 
dyprime.dx0) 

dx.dyO <- (cos.theta * dxprime.dy0) - (sin.theta * 
dyprime.dy0) 

dy.dyO <- (sin.theta * dxprime.dy0) + (cos.theta * 


dyprime.dy0O) + 1 


Grad. <0) <2 e Sum x-a1 tf ~fdx.dx0 
Gadel yO tee SU (x. adltte- ax. dy0 
. theta, 


Grad. <—— ,C\(Grad.a, 
grad.delta) 


Grad.b,) grad 


Ie) {| 


there.hess == 
print (grad) 


return (grad) 


Pied S 


} 

d2xprime.da2 <- d2yprime.da2 <- 0 
d2xprime.dade <- dxprime.de/a 
d2yprime.dade <- dyprime.de/a 


+ 
~- 


Y.diff£i* 'ayvdz0} 
Voadl1ff * sdyeay0} 
Grad. x07. Grad. 70, 


d2xprime.dadtheta <- dxprime.dtheta/a 
d2yprime.dadtheta <- dyprime.dtheta/a 


G@denom.de <- -2Z2 * e * cossg.tt.theta 
ddenom.dtheta <- - e*“2 * sin(2 * 
berm <— 
xprime.denom.sq <- xprime * denom*2 


d2xprime.de2 <- xprime.denom.sq - e * 


xprime * denom * ddenom.de) 
d2xprime.dez <- 
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(terml/xprime.denom.sq%2) 


(tt -— theta) ) 
(="a°Z)*. Sinsq=2-ct.theta)/4 


(axprime.der* denom 2 + 2 ~ 


* d2xprime.de2 


HH # 


HH # 


HH # 
HHH 


HHH 
# # # 


d2xprime.de2[abs(xprime) < le-006] <- 0 


d2yprime.de2 <—- one minus-e.sq * {(-2 * xprime *~ d2Zxprime.de2Z — 72a 
dxprime.de*2) + 8 * e * xprime * dxprime.de - 2 * (a%2 - 
xprime%2) 

terml <- - a*%2 * e 

d2xprime.dedtheta <- - xprime.denom.sgq * sin(4 * (tt - theta))} = 


Sinsgq.2.tt.theta * (xprime * denom * ddenom.dtheta + 
dxprime.dtheta * denom’*2) 

d2xprime.dedtheta <- (( - a*2 * e)/xprime.denom.sq%*2) * 
d2xprime.dedtheta 

d2xprime.dedtheta[abs(xprime) < le-006] <- 0 

d2yprime.dedtheta <- (-1/yprime*2) * (yprime * ((one.minus.e.sq * 


xprime * d2xprime.dedtheta + dxprime.dtheta * dxprime.de)) - 
2* e * xprime * dxprime.dtheta) - one.minus.e.sq * xprime * 
dxprime.dtheta * dxprime.de) 

d2yprime.dedtheta[abs(yprime) < le-006] <- 0 

dnum.dtheta <- -2 * one.minus.e.sgq * cos(2 * (tt - theta) ) 

d2xprime.dtheta2 <- xprime.denom.sgq * dnum.dtheta - num * (2 * 
xprime * denom * ddenom.dtheta + denom*’2 * dxprime.dtheta) 

d2xprime.dtheta2 <- (d2xprime.dtheta2 * a%*2)/(2 * 
xprime.denom.sq*%2) 

d2xprime.dtheta2[abs(xprime) < le-006] <- 0 


d2yprime.dtheta2 <- - (one.minus.e.sgq/yprime) * (yprime * (xprime 
* d2xprime.dtheta2 + dxprime.dtheta%*2) - dyprime.dtheta * ({ 
xprime * dxprime.dtheta) ) 

d2yprime.dtheta2[abs(yprime) < le-006] <- 0 4 


d2e.db2 <- (b * de.db - e)/(a * e)%2 

d2e.dadb <- (2 * b)/(a*3 * e) # 

d2x.dade <- cos.theta * dZ2xprime.dade - sin.theta * d2Zyprime.dade 

d2x.dadb <- cos.theta * (dxprime.de * d2e.dadb + d2xprime.dade * 
de.db) - sin.theta * (dyprime.de * d2e.dadb + d2yprime.dade 
* de.db) # 

d2y.dade <- sin.theta * d2xprime.dade + cos.theta * d2yprime.dade 

d2y.dadb <- sin.theta * (dxprime.de * d2e.dadb + d2xprime.dade * 
de.db) + cos.theta * (dyprime.de * d2Z2e.dadb + d2yprime.dade 
* de.db) # 

d2x.de2 <- cos.theta * d2xprime.de2 - sin.theta * d2yprime.de2 

dZy.dezZ <-— sin-theta * dZxprime.deZ2 +Weos.theta *“#dZyprame .deZz 

d2x.dadb <- cos.theta * (dxprime.de * d2e.dadb + d2xprime.dade * 
de.db) - sin.theta * (dyprime.de * d2e.dadb + d2yprime.dade 
* de.db) 

d2y.dadb <- sin.theta * (dxprime.de * d2Z2e.dadb + d2xprime.dade * 
de.db) + cos.theta * (dyprime.de * d2e.dadb + d2yprime.dade 
* de.db) 

d2xprime.dedb <- d2xprime.de2 * de.db 

d2yprime.dedb <- d2yprime.de2 * de.db 

d2x.db2 <- cos.theta * (dxprime.de * d2e.db2 + d2xprime.dedb * 
de.db) - sin.theta * (dyprime.de * d2e.db2 + d2yprime-dedoa 
de .db) 

d2y.db2 <- sin.theta * (dxprime.de * d2e.db2 + d2xprime.dedb * 
de.db) + cos.theta * (dyprime.de * d2e.db2 + d2yprime.dedb * 
de. db) 

Grad.a2 <-> 2 = Sumi(d+-cdaecn cy 4aa-2) # 

Grad.ae <= 2 * sum( = x.diti *~ "d2x dade ax .da 1ax.de"—-"y.d1i rues 
d2y.dade + dy.da * dy.de) 

gGrad-ab <——25* sum@— x-ditt *@e2x.dadb + dxada * cieweb =| Y.ditr ae 
dZy.dadb==dy.da *#ay db) 
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# # # 
it H # 


# HH 


# # # 


# tH 
# HH 


SE SE OE OE 


eaZe Gad erePana wcoOs.theta ~ dzZXprime.dadtheta - sinitheta * 
d2yprime.dadtheta - dy.da 

d2y.dadtheta <- sin.theta * d2xprime.dadtheta + cos.theta * 
d2yprime.dadtheta + dx.da F 

d2x.dedtheta <- cos.theta * dZxprime.dedtheta - sin.theta * 
dZ2yprime.dedtheta - dyede 

d2y.dedtheta <- sin.theta * d2xprime.dedtheta + cos.theta * 
d2yprime.dedtheta + dx.de 

d2x.dbdtheta <- cos.theta * d2xprime.dedtheta * de.db - sin.theta 
* d2yprime.dedtheta * de.db - dy.db 

d2y.dbdtheta <- sin.theta * d2xprime.dedtheta * de.db + cos.theta 
* d2yprime.dedtheta * de.db + dx.db 

d2x.dtheta2 <- cos.theta * d2xprime.dtheta2 - sin.theta * 
d2yprime.dtheta2 - 2 * dy.dtheta + x 

dzy.dthetaZ2 <— sim-eheta * dzZxprime.dtheta2 + cos:.theta * 
dZyorime-dtheta2 + 2 *"dx.dtheta + y 

G@rad-aeteca <— 2 esum( — X.diff * dZx.dadtheta + dx -da * 

Gx GQeEhnela smyedlit ~ ad2y dadtheta + dy.da * dy.dtheta) # 

Grad. c2< en SUmin— xX dlrr ee dzZx<ade2 + dx.de“2 — y.diff * ### 
dZy.deZ + dytde~2) 

grad.etheta <- 2 * sum( - x.diff * d2x.dedtheta + dx.de *dx.dtheta 

- y.difi * d2y.dedtheta + dy.de * dy.dtheta) 

Gisele x. Ohba Ax .db2°+ ax.db°2.— y.diff * 

day.db2 + dy.db~2Z) 

grad. Dtheta<—eew Sumt( — xX-d1ff ~“d2x dbdtheta + dx.db * 
dum@etheta — y.diftf * dZy.dbdtheta + dy.db * dy.dtheta) 

Gada etetazea. 2 A GoUM( gaa x ditt “~“sd2x dthetaZ + dx-dtheta*2 — 

y-.diff * d2y.dtheta2 + dy.dtheta%2) 

if(fit.center == F) 
hessian <- c(grad.a2, grad.ab, grad.b2, grad.atheta, 

grad.btheta, grad.theta2) 
else { 


Second derivatives: a and x0, a and y0 


d2xprime.dadx0 <- dxprime.dx0/a 

d2yprime.dadx0 <- dyprime.dx0/a 

a2xprime.dadyO <- dxprime.dy0O/a 

d2yprime.dadyO <- dyprime.dy0/a 

d2x.dadx0 <- cos.theta * d2xprime.dadx0 - sin.theta * 
d2yprime.dadx0 

d2y.dadx0 <- sin.theta * d2xprime.dadx0 + cos.theta * 
d2yprime.dadx0 7 


# and the gradient 


GradtaxOs— aa SUNVGr- xX Gli teres. dadx0=+ dx-da * dx-dx0i— 
Voit reed cdadx Otay sca ~ dy. dx) 

d2x.dady0O <- cos.theta * dZ2xprime.dadyO - sin.theta * 
d2yprime.dady0 

d2y.dadyO <- sin.theta * d2xprime.dady0O + cos.theta * 


d2yprime.dady0 # 
grad.ayO <- 2 * sum( - x.diff * d2x.dadyO + dx.da * dx.dy0o - 
y.ditfis* dzZyedadyo s+ dy.da ~*~ dy.dy0) m 


ddenom-dxG@™=<=—e°2 * sin. Z2-tt.theta’* dt. dxo 
ddenémsdvee<- eo2°* Sin.2.tt.theta ~ dt.dyo 
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# 


# # # 
+H # 
HHH 
HHH 
Ht 
HH # 


HH 
HHH 
tH # 
HHH 


# 


A <= xprime * denom”2 


anh dx0' <= denom = (2 = sxprime *“ddenom-dx0 + denom = 


dxprime.dx0) 


GA.dy0O <- denom * (2 * xprime * ddenom.dy0O + denom * 


dxprime.dy0) 


out.front <- - (a*2 * e)/4 


z<- Sinsq.2.tt.theta/A # 


z{[abs (A) < tol] <- 0 


d2xprime.dedx0O <~- out.front * 


dt.dx0)/A - ((dA.dx0O * z)/A)) 


d2xprime.dedyO <—- out.front * 


dt.dy0)/A - ((dA.dyO * z)/A)) 
d2xprime.dedxOf[abs (A) < tol] <- 0 
d2xprime.dedyO[abs(A) < tol] <- 0 # 
# Here's one from Mathematica. 


aZzvyerLme.aca<6) <——(s— 


sin.2.tt.theta) /denom%*2 
dzyprime.dedxO[/abs(yprime) < tolje<— © 
d2x.dedx0 <= cos.theta * d2xprime.dedx0 - sin.theta * 


d2 yprime .dedx0 


((2>*ssin(4 * (Gt = theta). 


((2 * -sin(4 * (tt — theta) ama 


(3° *%5e)/2.* yorime * dtsdx0 > 


# 


G2 y dedxGm—— sinetnetals GdZxprrme aeadx0) + cos. theta * 


d2yprime.dedx0 
Grad vexOe<— 92 * sum(a— 


x. 


Y-Gier * a2y .dedxo 


dZ=. abax0.-—scos,cheta 
* d2yprime.dedx0 
G2y> dbax0=-= sin. theta 
* d2yprime.dedx0 
grad.bx0 <- 2 * sum( — 


* 
* 
* 
* 


as 


y.diff * d2y.dbdx0 


d2yprime.dedyO <- ( - 


d2yprime.dedy0 
d2y.dedy0 <- sin.theta 
d2yprime.dedy0 
d2x.dbdy0O <- cos.theta 
* d2yprime.dedy0 
d2y.dbdy0O <- sin.theta 
* d2yprime.dedy0 
GradebyC <—s2. sunt 


y-diff * d2y.dbdy 


# Here's another from Mathematica. 


# 


out.front <- (xprime * 


d2xprime.dthetadx0 <- out.front * dt.dx0 
dZ2xprime.dthetadyO <- out.front * dt.dy0 


eS 


* 


* 
“* 
* 
~ 
x 
0 


aqiff * d2x.dedx0 + dx.de * dx.dx0O - 


Peay sae. ay.cdx0) 


d2xprime.dedx0 * de.db - sin.theta 


de.db 


d2xprame.dedx0) = dé.db + cos. encima 


de.db 


Gift * d2x.dbdx0 + dx.db * dx 2d 


+ dy.db * dy.dx0) 


“eye "= Verime ~den,dy0u. 
sin.2.tt.theta) /denom*2 

d2yprime.dedyO[abs(yprime) < tol] <- 0 

d2x.dedy0 <- cos.theta * d2xprime.dedy0 


d2xprime .dedy0 


d2xprime.dedy0 
ae.dab 
d2xprime.dedy0 
ade.ab 


+ 


* 


* 


.diff * d2x.dbdy0 
+ dy.db * dy.dy0) 


# 


Sin.theta * 
cos.theta * 
de.db - sin.theta 
de.db + cos.theta 


+ dx2db.* dx -dytie 
# 


(1 - 2 * e*%2 + e%2 * 
cos.2.tt.theta))/ denom%*2 
out.frontfabs(denom) < tol] <- 0 


out.front <= ((one.minus.e.6q) * yvprime 


denom*2 


out.front{abs(denom) < tol] <- 0 


d2yprime.dthetadx0 <- out.front * dt.dx0 
d2yprime.dthetady0O <- out.front * dt-.dy0 
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7G 


# 
(3 - 2 * denom))/ 


# 


+t # 


dee eeledetd UN GCyecxO + cos.theta * d2xprime.dthetadx0 
- Sin.theta * d2yprime.dthetadx0 
Sgyounectadnm<-s(ax.ax0 — 1) + sin.theta * 
d2xprime.dthetadx0 + cos.theta * d2yprime.dthetadx0 
Grad thetaae) <— 2a Sim — x.difies d2x.dthetadx0 + 
dx clenbetay “dx. dxO07= y.d2££ * d2y.dthétadx0 + 
Gyedtheta * dy.dx0) 
CZ2esGenecrolad/0.——mm—mmcyocdyO - 1) + cos.theta * 
d2xprime.dthetady0 - sin.theta * d2yprime.dthetady0 
i2yecenerady0Os-— ax ayo + sin.theta * d2xprime.dthetadyo + 
cos.theta * d2yprime.dthetady0 
SaeeneeavOucoee = — stim, — X-Ql£f * dzZx.dthetadyO + 
Seearnetam sox avO = y.diff = dzy.dthetadyo + 
dy.dtheta * dy.dy0) 
# 


Cee Oni (OE ax0..* dt.dy0) 

d2t.dxO0dy0O <- -1/R.sq + 2 * (dt.dx0)%*2 

Zt 402 9 —  — dz2t.dx0Z 

d2xprime.dx02 <- - dxprime.dtheta * d2t.dx02 - 
d2xprime.dthetadx0 * dt.dx0 

d2yprime.dx02 <- - dyprime.dtheta * d2t.dx02 - 
dZyprime.dthetadx0 * dt.dx0 

d2x.dx02 <- cos.theta * d2xprime.dx02 - sin.theta * 
d2yprime.dx02 

aZy-ax02Z <— sin-theta * dZxprime.dx0Z + cos.theta * 
d2yprime.dx02 

SrA Veeco eee SUM (e= XeGtEt * dZ2x.dk02 + dxvdx0°2 = sy edirtet 
niel2yecxe2Z at dyodx0 72) 


d2xprime.dxOdyO <- - dxprime.dtheta * d2t.dxOdy0O - 
d2xprime.dthetadyO * dt.dx0 

dZypErime- ax0dy0e<—-S=—dyorime.dtheta * dZt.dx0dy0 — 
d2yprime.dthetady0O * dt.dx0 

d2x.dxO0dy0O <- cos.theta * d2xprime.dx0dy0 - sin.theta * 
d2yprime.dx0dy0 

adzyaeaxzeayO <= Sinstheta *~ “d2xprime.dx0dy0 + Cos:.theta * 
d2yprime.dx0dy0 

Geade 00s <a esiim( — x diff * A2xk dxz0dy 0b ds ax0 = 
dxadyOe==yoditt * d2y-dx0dy0™= 7 dy sao “sdy-dy0) # 


d2xprime.dy02 <- - dxprime.dtheta * d2t.dy02 - 
d2xprime.dthetady0O * dt.dy0 
d2yprime.dy02 <- —- dyprime.dtheta * d2t.dy02 - 


d2yprime.dthetady0O * dt.dy0 

d2eeavy02 <= cos-thneta * dZxpraime.dy02 — sin.theta * 
d2yprime.dy02 

d2y.dy02 <- sin.theta * d2xprime.dy02 + cos.theta * 
d2yprime.dy02 

GiadayU2.c-2 ou =k Crh er eee ayO2 + GX.dV00 2 = Vy clet 
Sra2y.ay02 1 ay .ady0- 2} 

hessian <- c(grad.a2, grad.ab, grad.b2, grad.atheta, 
grad.btheta, grad.theta2, grad.ax0, grad.bx0O, 
Grad. Ehetax0, grad.x02, -grad-ay0, grad-byG, 
grad.thetay0, grad.x0y0, grad.y02) # 

hessian <- c(grad.a2Z, grad.ae, grad.e2, grad.atheta, 
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# # # G@rad-etheta, .grad-thetaZ, grad.ax0,.qrad exo, 


# # # Grad ehetax0, grad. x02, grade ayO ms Grau. evo, 
# # # grad.thetay0, grad.x0y0, grad.y02) # 
# # # print (hessian) 


} 
thing <- list(gradient = grad, hesSian = hessian) 
return ching) 


> ell.pred 

function(tt, a, b, theta = 0, return: unrotated.too = F, fit.center =em 
center.x = 0, center.y — 0) 

{ 

# 

# Get fitted x and y for ellipse with points at angles "tt", 

# with eccentricity = e and a = a. Rotate by theta afterwards, 

# if asked. Finally, if fit.center = T, move everything by 

# center.x in the x direction and by center.y in the y direction. 

# 

# A little algebra shows that 

# a“2 Sin@2.pw 2 — tej (1, - 642) 

se ee — eh tau 

# Sin” 2( tt) tesin-Z (pi cee tt) eel —.e- 7) 

# 

# (If a < b, that's y*2, except you have to switch the tt's and the 

# (pi/2 = tts ee ihe sin 2 (pi, 2-ce)etermeais  ‘thang..) coptake 

# x (if a > b) to be the positive square root of that for the 

# moment. Then y*2 = (a - ex)*%2 - (ae —- x)*%2. So get that, too. 

# 


new.tt <- tt - theta 
if(a Seb) 


e <- sgqrt(l - (b/a)%*2) # 

## # if(e > 0.99) return(1000 * length(x) ) 
thang <- (sin(pi/2 - new.tt)*2) * (1 - e%2) # 
K <- SGEE( (as2e" = chang) (Ssininew.tt) 2 + thang).) 
yy <- (a - e * x)*2 - (a * e - x)*2 # 


# Make sure y*2 is always positive (round-off errors can hurt here); 
# then get y. 


ij 
VY [Y¥eeee | <— 50 
Yo<= SGruiyy) y 
} 
else { 
e-<— sgrt(l - (a/bje2) 4 
HF if(e > 0.99) return(1000 * length(x) ) 
thang <- (Sim(newete) -2Zie* Vibe— eZ) # 
yY <— by sqrte(thang/ (san (pi7 27] new- tt) 2+ Chang) ) 
xX <- (b - e * y)*%2 - (b * e - y)%2 # 
Sexxy Ol <—. oO 
x <- sgrt (xx) = 


) 
quad <=" newote.ss a2. pi) 
quad.2.3 <= quad > 2/2 & GQuads<cuo ply Z 


x[quad.2.3] <- - x[quad.2.3] 

quad . 3.74 sauce apa 

y (quad. 3-4) "<—" =) vilemad= 324) 

rotated.data <- matrix(c(cos(theta), - sin(theta), sin(theta), 
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COs\CMetNcCtai),, Za ee os rbind(x, y) 
1f(fit.center == F) { 

center.x <- 0 

center.y <- 0 


} 


TE{(KCEUrN.UnEocated.too == T) 
return(list(x = rotated.data[l1, | *-centersxyey = 
retated.datafZ, ] + center.y, x.prime = x, y.prime 
es) 
else VEclrumimgLLsin| Xe=eskotated datal[?t, )] + center.x, y = 


rotated.data([2, fer center.y)) 
} 


eeeLl .tt 

menction(x, y) 

{ 

# ell.tt: get angle values for data. It's acos(x/r) 

# in the first quadrant, etc. 

# 
tt <- numeric(length (x) ) 
Pat Lom<— x7 sqrt (xo2 + y*2) 
Baclo(ratio.> 1) -<— 1 
ratio[ratio < -1] <- -l 
ind <- x >= 0 & y >= 0 
tering) <= acos(ratio(ind] } 
ind <- x < 0 & y >= 0 


tt[ind] <- pi - acos( - ratio[ind]) 

ind <- x <0 &y < 0 

Eertnal <= pi + acos( = ratio(lind]) 

ind <- x >=0O0 & y < 0 

PeLIndi 9 <— 2 *~ pil = acos(ratzol[ind)) # 
EL 


Gy 





APPENDIX B. EXAMPLES 


pelssl = pc is subject 1’s initials 
the 1 is the spatial frequency 
ss1 means sum of squares obtained from ellipse (vice ellipse. or 
ellipse. 1) 

pelall = a data frame containing all of the data at 1 cpd for pc 
this data always consists of non-oblique data followed by oblique data. 
For this subject the first 80 (x,y) pairs are non-oblique and 81-160 are 
oblique data 

fit.center=T lets the ellipse center “float” vice being pinned to the origin 

> pclssl_ellipse(pclall[,1], pclall[,2], fit.center=T) 

a: 0.0615 ,b: 0.05194 ,th: -0.005236 ;x,y: 0.001551 0.00008819 ;obj: 0.06455 

a: 0.02742 ,b: 0.02742 ,th: -0.07295 ;x,y: 0.01894 -0.001528 ;obj: 0.07366 

removed middle output to save space 

a: 0.04698 ,b: 0.03466 ,th: -0.02604 ;x,y: 0.004495 0.0006379 ;obj: 0.02068 

a: 0.04698 ,b: 0.03466 ,th: -0.02605 ;x,y: 0.004495 0.0006379 ;obj: 0.02068 

a: 0.04698 ,b: 0.03466 ,th: -0.02605 ;x,y: 0.004495 0.0006379 ;obj: 0.02068 

> pclssl 

Message was 

RELATIVE FUNCTION CONVERGENCE 
a b ~—stheta’=«scenter.x _—scenter.y 

0.046983 78 0.034663 18 -0.02604936 0.004495335 0.000637888 1 

> pclssl[2] 

Sobjective: 

[1] 0.02068425 


pcelss2. 01 = pc is subject 1’s initials 
the 1 is the spatial frequency 
ss2 means sum of squares obtained from ellipse. or ellipse. LI 
-III lets us know this is from ellipse. 
> pclss2.IIL_ ellipse.III (pclall[,1], pclall[,2], grad=F, 
+ is.there.hess=F, fit.center=T, class.I = (1:160) < 81, plot=F) 
a: 0.07282 , delta.a: 0 ,b: 0.0615 delta.b: 0 ,th: -0.005236 ;x,y: 0.001551 
0.000088 19 ;obj: 0.1364 
a: 0.07283 , delta.a: 0 ,b: 0.0615 delta.b: 0 ,th: -0.005236 ;x,y: 0.001551 
0.00008819 ;obj: 0.1364 
removed middle output to save space 
6 ;x,y: 0.004467 0.0006382 ;obj: 0.02067 
a: 0.04703 , delta.a: -0.0001285 ,b: 0.03424 delta.b: 0.0008527 ,th: -0.0266 
6 :x,y: 0.004467 0.0006382 ;obj: 0.02067 


71 


a: 0.04703 , delta.a: -0.0001281 ,b: 0.03424 delta.b: 0.000853 ,th: -0.02666 
‘x y: 0.004467 0.0006382 ;obj: 0.02067 
a: 0.04703 , delta.a: -0.0001281 ,b: 0.03424 delta.b: 0.0008523 ,th: -0.0266 
6 :x,y: 0.004467 0.0006382 ;obj: 0.02067 
Warning messages: 
1: singularity encountered in: nlminb.0(temp, p, liv, lv, objective, bounds, 
scale) 
removed identical warnings numbered 2 and 3 
4: singularity encountered in: nlminb.O(temp, p, liv, lv, objective, bounds, 
scale) 
> pclss2.II 
Message was 
RELATIVE FUNCTION CONVERGENCE 


a b theta center.x  center.y delta.a 
0.04703 421 0.03424435 -0.0266577 0.004466774 0.0006382473 -0.0001281251 
delta.b 
0.0008526696 


this is the objective function value which the program minimizes 
> pclss2.III[2] 

Sobjective: 

[1] 0.02067104 


p-values for the ellipses being different? 

> 1-pf(((pelss1 [[2]]-pelss2.T0[[2]])/2)/(pelss1[[2]]/(160-5)),2,155) 
[1] 0.9517269 

> 1-pf(((pce3ss1[[2]]-pc3ss2.T0[[2]])/2)/(pce3ss1[[2]]/(1 60-5)),2,155) 
[1] 0.07660905 

> 1-pf(((pe7ss1[[2]]-pc7ss2.TI[[2]])/2)/(pc7ss1[[2]]/(160-5)),2, 155) 
[1] 1.822829e-007 


are the b’s significant 

> 1-pf((pelss2.II[[2]]-pcelss2.1II[[2]])/(pclss2.1M[[2]]/153), 1,153) 
[1] 0.7693711 

> |-pf((pce3ss2.II[[2]]-pce3 ss2. III [[2]])/(pe3ss2 IN[[2]]/153), 1,153) 
[1] 0.837154 

> 1-pf((pe7ss2.II[[2]]-pc7ss2.1II[[2]])/(pe7ss2.M1[[2]]/153), 1,153) 
[1] 2.110824e-008 


changing which.type to 2 forces the non-oblique and oblique ellipses to have the 
Same major axis. 


da 


> pclss2.I1.which.type.2_ ellipse. II(pclall[, 1],pclall[,2],fit.center=T,class.I = (1:160) < 
8 1,which.type=2,grad=F,1s.there.hess=F , plot.it=F) 

>pc3ss2.II.which.type.2_ ellipse. I(pc3all[, 1],pc3all[,2],fit.center=T,class.I = (1:160) < 

8 1,which.type=2, grad=F ,1s.there. hess=F, plot.it=F) 

>pc7ss2.I].which.type.2_ ellipse. I(pc7all[, 1],pc7all[,2], fit.center=T,class.I = (1:160) < 

8 1,which.type=2,grad=F,1s.there.hess=F, plot.it=F) 


are the a’s significant ? 

subject 1 at 7 cpd 

> 1-pf((pc7ss2.II.which.type.2[[2]]-pc7ss2.III[[2]])/(pce7ss2. malay 13) lid 3) 
[1] 0.5883988 

subject 1 at 3 cpd 

> -pf((pce3ss2.II.which.type.2[[2]]-pc3ss2.II[[2]])/(pce3ss2.III[[2]]/153), 1,153) 
[1] 0.02835718 

subject 1 at 1 cpd 

> 1-pf((pelss2.I..which.type.2[[2]]-pclss2.U1[[2]])/(pcelss2.IM[[2]]/153), 1,153) 
[1] 0.9753455 
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